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1.  Introduction 

1.1.  Background 

Process  and  device  design  for  modem  integrated  circuits  depends  on  a  complex  interrelationship 
of  manufacturing  techniques,  physical  phenomena,  device  requirements,  and  designer  experience.  With 
increases  in  integrated  circuit  complexity  and  decreases  in  device  feature  dimensions,  computer-aided- 
design  software  has  come  to  play  an  important  role  in  product  development.  It  is  now  practically  incon¬ 
ceivable  that  an  integrated  circuit  could  be  designed  without  the  help  of  circuit  simulators  such  as 
SPICE  [Nag75]  and  layout  routing  tools  such  as  TimberWolf  [Sec85].  Most  commercial  enterprises 
involved  in  chip-making  have  robust  software  and  computing  environments  available  for  IC  simulation, 
routing,  and  design.  Such  is  however  not  the  case  for  process  and  device  modeling.  Despite  the  many 
successes  of  device  simulators  such  as  PISCES  [Pin84]]  and  MINIMOS  [Sel80],  and  process  simulators 
such  as  SAMPLE  [01d79]  and  SUPREM  [Ho83a],  process  and  device  simulation  is  not  nearly  as  well 
established  in  industry  as  circuit  simulation  and  layout  automation.  This  is  evidenced  by  the  fact,  that 
whereas  it  is  now  nearly  impossible  to  design  a  circuit  without  a  circuit  simulator,  it  is  still  possible  for 
device  and  process  engineers  to  do  their  jobs  without  availing  themselves  of  the  software  capabilities  at 
hand. 

Process  and  device  modeling  is  a  difficult  task  given  the  variety  and  complexity  of  physical 
phenomena  which  must  be  considered.  Even  when  the  physics  of  a  certain  situation  can  be  reduced  to  a 
handful  of  equations,  often  algorithms  for  numerical  solution  require  computing  power  beyond  that 
readily  available.  The  challenge  for  process  and  device  simulation  involves  developing  accurate 
models,  determining  parameter  values  for  those  models,  and  then  implementing  algorithms  which  per¬ 
form  the  necessary  calculations  in  a  reasonable  amount  of  time.  To  achieve  these  goals,  many  existing 
simulators  depend  on  simplifying  assumptions  which  restrict  their  applicability  and  accuracy. 
SUPREM-Hl  [Ho83b]  only  handles  diffusion  and  oxidation  in  one  dimension.  SAMPLE- 1.6a  [Add85] 
requires  planar  layers  for  lithography  and  etching  simulation  in  two  dimensions.  Additionally,  many  of 


the  physical  models  used  by  process  simulators  are  heavily  dependent  on  empirical  results.  SAMPLE, 
for  example,  requires  experimentally  derived  development  rate  equations  and  curve-fitting  parameters 
for  each  photoresist  material  it  is  required  to  simulate.  The  empirical  nature  of  many  of  the  models 
limits  the  range  over  which  input  parameters  may  be  varied  and  thus  limits  the  usefulness  of  the  simu¬ 
lation  software. 

New  versions  of  these  programs,  however,  are  closing  the  gap  between  process  simulation  and 
process  reality.  SAMPLE-1. 7a  [SAM88]  includes  models  for  non-planar  etching,  and  advanced  resist 
chemistries.  A  program  for  development  rate  parameter  extraction,  PARMEX  [Bel88],  will  soon  be 
available  for  rapidly  characterizing  photoresists  for  SAMPLE.  SUPREM-IV  [Rafi86]  handles  diffusion 
in  two  dimensions.  CREEP  [SutS7]  is  a  new  program  for  oxidation,  annealing  and  reflow  simulation. 
COMPOSITE  [Lor85]  includes  many  improved  etching  models  as  well  as  providing  an  extensive  library 
of  topography,  diffusion  and  implantation  simulation  programs. 

Advances  in  process  simulation  are  moving  in  the  direction  of  complete  three-dimensional  model¬ 
ing.  Already,  the  program  SPLAT  [Toh88]  is  capable  of  simulating  the  variation  of  light  intensity 
throughout  a  field  on  a  wafer  in  to  a  microlithography  system,  including  the  effects  of  lens  aberrations. 
Future  work  will  involve  simulating  resist  dissolution  and  pattern  transfer  in  three  dimensions  in  order 
to  give  a  more  complete  physical  picture  of  process  phenomena.  The  advances  in  physical  modeling 
for  process  simulators,  coupled  with  the  recent  availability  of  powerful  "super-minicomputers"  indicates 
that  this  trend  towards  improved  process  simulation  will  continue. 

1.2.  Problem  of  Process  CAD  Integration 

One  problem  that  has  arisen  as  a  direct  result  of  the  increased  number  of  simulation  programs 
available,  is  the  difficulty  of  passing  data  from  one  simulator  to  another.  Production  line  processes 
depend  on  the  interaction  of  several  different  processing  steps.  Existing  process  simulators  generally 
concentrate  on  only  a  certain  fraction  of  possible  processes.  For  many  research  and  design  applications 
this  is  sufficient,  but  there  is  now  a  more  pressing  need  for  integrating  the  many  available  CAD  pro¬ 
grams  for  process  simulation.  It  is  also  desirable  to  use  process  simulators  to  create  input  data  for 


device  simulators  so  that  the  electrical  characteristics  of  a  given  device  can  be  examined  with  respect  to 
process  variations.  Additionally,  process  simulators  should  be  connected  to  the  layout  so  that  the  effect 
of  design  rules  and  mask  alignment  on  topography  can  be  considered. 

The  latter  problem  of  connecting  the  layout  to  the  simulator  has  already  been  addressed  by 
SIMPL-1  and  SIMPL-2  [Gri84],  [Lee85].  Both  of  these  programs  were  designed  to  generate  a  device 
cross  section  along  a  cut-line  on  the  layout  based  on  a  process  flow  description.  SIMPL-2  can  represent 
a  device  cross  section  using  arbitrary  polygons  and  is  therefore  quite  general  when  it  comes  to  describ¬ 
ing  a  device  geometry.  SIMPL-2  lacks  sophisticated  process  models,  however,  and  is  thus  of  primary 
interest  to  layout  designers  but  not  to  process  engineers.  SIMPL-2  offers  rigorous  simulation  only  in  the 
case  of  metal  deposition,  in  which  case  SAMPLE  is  called.  SIMPL-2  also  lacks  many  convenient 
features  from  a  user’s  point  of  view.  The  graphics  are  slow  and  it  is  not  possible  to  manipulate  the 
display  easily. 

SIMPL-DIX  [Wu88a]  solves  many  of  the  problems  of  graphics  and  user  friendliness  by  providing 
a  design  interface  using  the  X  window  system  [Sch86].  SIMPL-DIX  also  supports  a  variety  of  features 
for  analyzing  the  effects  of  mask  misalignment  (WORST),  problems  due  to  the  layout  and  geometric 
effects  (HUNCH),  and  characteristics  of  the  simulated  device  profile  itself  (CRITIC).  Furthermore, 
SIMPL-DIX  combined  with  SIMPL-2  is  well  positioned  to  fill  the  need  for  an  integrated  design 
environment,  providing  for  the  flow  of  information  between  process  simulators  as  well  as  acting  as  a 
front-end  input  generator  for  device  simulators.  Some  initial  work  in  integrating  the  SIMPL  programs 
with  CREEP,  SUPREM-III  and  PISCES-II  has  recently  been  reported  [Wu88b].  Figure  1.1  presents  a 
summary  of  existing  and  proposed  tools  integration  using  SLMPL. 

1J.  Project  Overview 

It  is  the  intent  of  this  project  to  demonstrate  the  feasibility  and  flexibility  of  SIMPL  as  an  inter¬ 
face  to  process  and  device  simulators.  The  usefulness  of  this  approach  will  be  demonstrated  as  SIMPL 
is  applied  to  a  few  selected  problems  suited  for  analysis  with  these  programs.  SIMPL  has  been  con¬ 
nected  with  RACPLE  [Lee83]  for  analyzing  topography  dependent  resistive  and  capacitive  parasitics. 
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SIMPL  is  also  used  for  displaying  the  results  of  non-planar  etch  simulation  using  SAMPLE- 1.7a.  This 
report  details  the  new  algorithms  and  interface  code  needed  to  integrate  RACPLE  and  SAMPLE  non- 
planar  etch  with  SIMPL.  These  integrated  CAD  tools  are  applied  to  problems  of  multilevel  planariza¬ 
tion  schemes  and  dense  interconnect  technologies.  Finally,  based  on  experience  derived  from  working 
on  this  project,  several  suggestions  for  improving  the  SIMPL  programs  are  offered. 
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2.  A  General  Approach  to  Software  Integration  with  SIMPL 

2.1.  SIMPL-DIX  and  SIMPL-2 

SIMPL-DIX  and  SIMPL-2  are  two  separate  programs  which,  together  with  SIMPL- 1,  form  a  suite 
of  programs  known  genetically  as  SIMPL  (SIMulated  Profiles  from  the  Layout).  It  is  worthwhile  to 
understand  the  capabilities  and  purposes  of  these  programs  in  a  general  way  before  considering  their 
specific  implementations.  SIMPL-2  was  originally  developed  to  allow  a  layout  designer  to  rapidly  deter¬ 
mine  the  cross  section  of  a  device  which  would  result  from  a  particular  layout.  To  achieve  this, 
SIMPL-2  uses  arbitrary  polygons  to  represent  device  cross  sections.  Crude,  but  fast  models  for  simulat¬ 
ing  deposition,  etching,  exposure,  development,  oxidation,  and  ion  implantation  are  available.  SIMPL-2 
is  capable  of  displaying  two-dimensional  process  effects  such  as  "bird’s  beak"  oxidation,  lateral 
diffusion  under  a  mask  edge,  and  undercut  in  etching.  A  link  to  SAMPLE  for  rigorous  simulation  of 
metal  deposition  is  included  to  describe  sidewall  and  step  coverage  accurately.  Using  SIMPL-2,  it  is 
possible  to  generate  a  realistic  cross  section  of  a  device  such  as  the  Berkeley  CMOS  inverter  [Gib86]. 

SIMPL-DIX  was  developed  to  provide  a  high-level  graphics  interface  and  convenient  user 
environment  for  programs  like  SIMPL-2.  A  major  goal  of  SIMPL-DIX  is  to  provide  an  environment  for 
integrating  dissimilar  process  and  device  simulation  programs.  SIMPL-DIX  uses  the  X  window  system 
for  displaying  graphics  information  and  for  generating  a  menu-driven  user  interface.  SIMPL-DIX 
currently  invokes  SIMPL-2  to  create  the  database  needed  to  describe  a  device  cross  section.  SIMPL- 
DIX  also  maintains  its  own  database  for  describing  all  polygons  in  a  profile,  but  this  database  is  pri¬ 
marily  intended  for  storing  graphics  information  and  thus  lacks  many  of  the  features  of  the  SIMPL-2 
database.  For  example,  SIMPL-DIX  does  not  maintain  data  describing  polygons  adjacent  to  a  given 
polygon,  whereas  SIMPL-2  contains  the  material  type  and  relative  location  of  all  polygons  sharing  a 
boundary  or  vertex  point. 

In  addition  to  the  link  with  SIMPL-2,  SIMPL-DIX  has  a  number  of  internal  tools  to  assist  the 
designer  in  performing  and  analyzing  simulations.  This  points  to  the  ability  of  SIMPL-DIX  to  analyze 
a  profile  or  layout  generated  by  other  programs.  The  HUNCH  feature  allows  a  designer  to  specify 
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logical  operations  between  masks  or  sets  of  mask  to  identify  locations  where  topographical  problems 
are  anticipated.  The  CRITIC  feature  is  currently  being  developed  to  allow  aspects  of  a  device  profile  to 
be  investigated  automatically. 

The  SIMPL  programs  together  form  a  suitable  design  environment  for  process  and  device  simula¬ 
tion.  SIMPL-DIX  provides  a  convenient  graphics  and  user  interface,  whereas  SIMPL-2  maintains  the 
most  complete  data  base  for  describing  device  profiles.  The  original  goal  of  SIMPL  to  link  process 
simulation  with  the  device  layout  can  now  be  extended  to  include  integration  of  other  device  and  pro¬ 
cess  simulators. 

22.  Profile  Data  Representations 

Several  schemes  exist  for  representing  device  cross  sections.  Most  can  be  classified  as  belonging 
to  either  one  of  two  types:  linked  polygons  or  multiple  layers.  SIMPL-2  uses  linked  polygons  to 
describe  a  device  profile.  This  technique  offers  the  advantage  of  being  quite  general.  Nearly  all  device 
profiles  can  be  conveniently  described  in  this  way.  Polygons  can  be  added,  removed,  or  altered  by  using 
a  suite  of  subroutines  to  manipulate  the  database.  Figure  2.1  shows  the  basic  setup  of  the  linked 
polygon  data  structure. 

An  alternative  approach  to  representing  device  cross  sections  is  to  use  multiple  layers.  This 
method  has  certain  computational  advantages  for  simulating  processes  such  as  etching,  resist  dissolu¬ 
tion,  and  deposition.  SAMPLE  and  COMPOSITE  use  multiple  layers  to  represent  device  profiles.  Gen¬ 
erally,  a  layer  must  span  the  entire  length  of  a  simulation  window.  In  regions  where  no  material  of  a 
given  layer  type  exists,  it  is  possible  to  specify  a  layer  of  zero  thickness.  Some  programs,  like  COMPO¬ 
SITE,  do  not  require  that  a  layer  span  the  entire  simulation  window,  and  instead  maintain  an  additional 
flag  at  each  vertex  in  the  layer  to  describe  whether  the  point  is  the  starting  or  ending  point  of  a  layer. 
Figure  2.2  shows  how  layers  can  be  used  to  describe  a  profile  in  SAMPLE. 

In  order  to  integrate  programs,  it  is  necessary  to  translate  data  between  different  data  representa¬ 
tions.  It  is  not  difficult  to  translate  a  polygonal  data  representation  into  a  series  of  layers.  This  can  be 
achieved  by  writing  the  surface  of  the  profile  as  one  layer,  removing  one  surface  polygon,  writing  the 


7 


new  surface  as  the  next  layer,  and  so  on  down  to  the  substrate.  Going  from  layers  to  polygons,  how¬ 
ever,  is  somewhat  more  difficult  Two  contours,  describing  one  polygon,  must  be  identified  and  merged, 
with  extraneous  points  being  removed.  Fortunately,  there  are  many  special  cases  where  it  is  not  neces¬ 
sary  to  translate  an  entire  profile  back  and  forth  between  different  simulators.  In  the  case  of  layer  depo¬ 
sition,  only  the  surface  layer  need  be  sent,  and  only  the  new  deposited  layer  must  be  returned.  In  the 
case  of  etch  simulation,  usually  the  entire  cross  section  must  be  sent,  but  only  the  new  resulting  surface 
contour  need  be  returned.  In  some  special  cases,  it  is  a  trivial  matter  to  include  this  new  contour  in  the 
original  profile.  In  other  more  general  cases,  all  polygons  above  the  new  contour  must  be  removed,  and 
polygons  that  intersect  it  must  be  clipped.  This  case  is  not  trivial  to  implement  Integrating  programs 
which  share  similar  databases,  such  as  COMPOSITE  and  SAMPLE,  or  SIMPL  and  CREEP,  is  concep¬ 
tually  simple  but  translating  among  profile  representation  formats  can  become  tedious  and  time- 
consuming. 

A  standard  Profile  Interchange  Format  (PIF)  has  been  proposed  [Duv88]  but  has  not  yet  been 
fully  implemented.  PEF  is  a  polygon  based  data  representation  which  could  be  used  for  passing  data 
among  different  programs.  It  would  still  be  necessary  to  translate  a  PDF  profile  representation  into  that 
used  by  a  particular  program,  but  the  difficulties  of  coordinating  different  formats  would  be  streamlined 
with  such  a  standard.  The  problem  of  tools  integration  would  be  greatly  eased  by  a  standard  format  for 
describing  device  cross  sections.  The  basic  concepts  of  data  translation  would  still  be  necessary',  but 
their  implementation  would  be  greatly  simplified. 

23.  SIMPL  as  a  User  Interface 

Another  area  to  be  addressed  involves  the  convenient  input  of  data  to  a  simulation  program. 
SIMPL  is  a  menu  driven  environment  which  greatly  facilitates  its  use.  A  user  need  not  be  intimately 
familiar  with  the  specifics  of  SIMPL  in  order  to  get  started  with  it  This  menu  feature  can  be  extended 
to  other  simulators  called  by  SIMPL.  Input  files  for  a  program  such  as  SAMPLE  can  be  built  up 
automatically  through  a  series  of  questions  presented  by  SIMPL.  The  less  time  someone  has  to  spend 
learning  how  to  use  a  program,  the  more  likely  that  person  is  to  use  it  If  less  time  is  spent  writing 
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input  files,  more  time  can  be  devoted  to  using  the  program  or  doing  other  work.  The  advantages  of  a 
menu  driven  user  interface  are  clear. 

SIMPL-DIX  offers  several  routines  for  creating  menu  options  on  the  display,  and  for  prompting 
data  input  from  a  user.  It  is  a  relatively  easy  task  to  create  new  menu  options  for  SIMPL  and  instruc¬ 
tions  for  doing  so  are  clearly  spelled  out  in  the  thesis  on  SIMPL-DIX  [Wu88a]. 

2.4.  Displaying  Data  with  SIMPL 

It  is  similarly  straight-forward  to  display  data  using  SIMPL.  All  of  the  data  structures  required  by 
the  X  window  system  are  already  set  up  in  SIMPL-DIX.  Routines  for  locating  points  in  the  display 
viewports  are  easily  used.  Basic  X  window  function  calls  can  be  used  to  display  numerical  data  on  the 
screen,  highlight  certain  regions  of  the  display,  or  draw  new  information  to  the  screen.  The  only  chal¬ 
lenge  in  displaying  data  with  SIMPL  is  deciding  precisely  where  on  the  screen  to  put  the  data  so  that  it 
is  most  useful  and  most  cleanly  represented.  Figure  2.3  is  an  example  of  displaying  data  associated  with 
individual  polygons  in  a  SIMPL  cross  section. 

2.5.  Unix  Interface 

SIMPL-DIX  uses  the  interprocess  communication  (EPC)  facilities  in  the  Berkeley  UNIX  4.3BSD 
release  [Sec85]  to  provide  a  connection  to  SIMPL-2.  This  method  of  communication  is  useful  for  pass¬ 
ing  information  to  an  interactive  program.  It  is  also  possible  to  communicate  with  programs  through 
files,  by  creating  an  input  file  for  a  given  program  and  then  using  a  Unix  system  call  to  run  that  pro¬ 
gram.  The  method  using  IPC  is  more  sophisticated  in  that  it  allows  multiple  process  to  be  run  simul¬ 
taneously.  Unix  system  calls,  however,  are  easier  to  implement  and  are  adequate  in  situations  where  it 
is  not  necessary  to  run  multiple  processes.  SIMPL  uses  both  methods  of  inter-program  communication. 
EPC  is  used  for  the  SIMPL-DIX  to  SIMPL-2  interface  since  both  programs  are  run  simultaneously  and 
short  streams  of  data,  such  as  process  commands,  are  continuously  sent  between  them.  Unix  system 
calls  are  used  to  run  SAMPLE  since  large  streams  of  data,  such  as  profile  contours,  must  be  sent  back 
and  forth.  Additionally,  nothing  is  gained  by  using  IPC  in  this  case,  since  SIMPL  must  wait  for  the  data 


from  SAMPLE  before  continuing. 
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3.  Interface  to  RACPLE  for  Parasitic  Evaluation  and  SAMPLE  1.7  for  Nonplanar  Etching. 

3.1.  The  RACPLE  Program 

RACPLE  (Resistance  and  Capacitance  of  Profiles  in  Lithography  and  Etching)  is  a  post  processor 
for  the  SAMPLE  program.  It  calculates  the  effective  ratio  of  length  vs.  depth  for  a  thin  film  deposited 
over  a  nonplanar  surface.  The  effective  number  of  lateral  squares  (length/depth)  can  be  used  to  deter¬ 
mine  the  resistance  along  the  film  or  the  capacitance  across  the  layer.  For  a  planar  conducting  layer,  the 
resistance  is  given  by 

R  =(resistivity  )(L  ID  )(1  IW ) 
the  capacitance  for  a  planar  dielectric  is  given  by 

C  =(j>ermittivity  )(L  ID)W 

RACPLE  calculates  the  effective  LID  for  a  thin  film  by  locating  critical  features  in  the  profile  and  then 
dividing  the  film  into  small  sections  for  which  the  resistance  can  be  approximated.  This  approach  has 
an  accuracy  of  better  than  5%  when  compared  with  a  numerical  solution  of  the  Laplace  equation.  Care 
should  be  taken  though,  since  RACPLE  does  not  include  the  effects  of  fringing  fields  for  determining 
capacitance.  RACPLE  should  only  be  used  for  capacitance  simulation  if  the  dielectric  is  thin  compared 
to  its  length  and  width.  Resistance  calculation  is  accurate  for  a  wide  range  of  geometries.  More  infor¬ 
mation  on  RACPLE  is  contained  in  the  original  report  on  the  program  [Lee83b]. 

RACPLE  is  consistent  with  the  basic  philosophy  of  SIMPL  in  providing  rapid  evaluation  of  a 
device  profile  with  a  minimum  of  computational  effort.  Despite  its  limitations,  RACPLE  is  a  useful  pro¬ 
gram  for  topography  dependent  electrical  parameter  extraction,  and  is  a  useful  addition  to  the  SIMPL 
design  environment. 

3.2.  SIMPL  -  RACPLE  Interface. 

SIMPL-DIX  is  used  to  generate  input  files  for  RACPLE,  call  RACPLE  to  calculate  parasitics,  and 
then  display  the  results  along  with  the  device  profile.  The  interface  extracts  each  polygon  in  the 
SIMPL-DIX  data  base,  creates  an  input  file  for  each  polygon,  and  then  runs  RACPLE. 
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RACPLE  uses  the  SAMPLE  plotting  file  format  (f77punch7)  as  its  input  data  representation.  The 
interface  takes  the  polygon  as  given  by  SIMPL-DIX  and  writes  it  as  a  top  contour  and  a  bottom  contour 
in  a  SAMPLE  format  plot  file.  The  polygon  as  represented  in  SIMPL-DIX  is  a  linked  list  of  coordinate 
pairs  which  traverse  the  polygon  boundary  in  a  clockwise  direction.  There  is,  however,  no  preferred 
location  for  the  first  vertex  in  the  list  The  entire  polygon  must  be  read  in  order  to  determine  which 
points  make  up  the  top  contour  and  which  make  up  the  bottom  contour.  Figure  3.1  gives  the  shape  of  a 
typical  generic  polygon.  The  left  and  right  edges  need  not  be  vertical,  but  this  is  often  the  case  for  a 
polygon  which  touches  the  edge  of  the  simulation  window.  The  interface  traverses  the  polygon  to  deter¬ 
mine  the  relative  location  of  the  first  vertex  in  the  list,  as  well  as  the  character  of  the  left  and  right 
edges.  The  head  vertex  is  classified  as  belonging  to  one  of  eight  types  as  shown  in  Figure  3.2,  either  a 
comer  point  0,  2,  4,  or  6  or  an  intermediate  point  1,  3,  5,  or  7.  From  this  information,  the  interface 
arranges  the  points  into  a  top  and  bottom  layer  which  are  written  to  a  file  in  SAMPLE  plot  file  format. 
Figure  3.3  shows  how  the  example  polygon  of  Figure  3.2  is  split  up  into  two  layers. 

The  routines  used  by  RACPLE  expect  a  high  density  of  points  to  describe  a  profile  contour,  even 
if  the  profile  is  a  line  segment.  SIMPL,  in  order  to  save  space,  eliminates  redundant  points  from  a 
polygon,  so  that  a  line  segment  is  described  only  by  its  end  points.  To  ensure  computational  accuracy 
with  RACPLE,  the  interface  inserts  additional  points  in  the  contours  until  no  pair  of  adjacent  points  are 
separated  by  more  than  about  l/20th  the  length  of  the  contour.  Failing  to  do  this  can  result  in  errors  of 
as  much  as  75%  for  typical  profiles,  since  RACPLE  cannot  split  up  the  profile  correctly  if  the  layers 
lack  a  sufficient  number  of  points. 

As  each  input  file  is  created,  it  is  run  with  RACPLE  and  the  results  added  to  a  file  called, 
appropriately  enough,  RACPLE_RESULTS.  The  interface  uses  the  Unix  system  command  to  call 
RACPLE  since  this  method  is  convenient  to  implement  and  sufficient  for  the  purposes  of  this  interface. 
Once  all  the  polygons  have  been  evaluated  by  RACPLE,  the  results  can  be  displayed.  The  interface  first 
determines  the  dimensions  of  the  display  viewport  using  global  variables  in  SIMPL-DIX.  Each  piece  of 
data  generated  by  RACPLE  is  associated  with  a  particular  polygon  which  has  one  point  in  its  upper  left 
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comer.  SIMPL-DIX  displays  the  RACPLE  data  along  the  top  and  bottom  of  the  profile  viewport  and 
draws  a  line  from  the  data  to  the  corresponding  polygon.  To  avoid  a  sloppy  display  with  many  crossing 
lines,  the  data  is  first  sorted  by  the  x-value  of  the  upper  left  vertex  of  the  polygon  associated  with  that 
piece  of  data.  For  polygons  that  have  upper  left  comers  with  the  same  x-coordinate,  an  additional  sort 
by.  the  y-coordinate  is  performed.  The  code  for  most  of  the  routines  used  in  the  RACPLE  interface  is 
listed  in  Appendix  B  as  file  dix_action4.c.  The  contents  of  this  file  are  compiled  as  part  of  SIMPL- 
DIX. 

RACPLE  is  called  from  the  CRITIC  menu  in  SIMPL-DIX  which  is  in  turn  part  of  the  TOOLS 
option  of  the  main  SIMPL-DIX  menu,  of  the  version  of  SIMPL-DIX  currently  in  use  in  our  research 
group.  RACPLE  was  included  as  a  CRITIC  option  because  it  is  a  program  for  criticizing  and  analyzing 
profiles,  and  thus  fell  naturally  into  that  category  of  programs. 

33.  SIMPL  interface  to  run  SAMPLE  1.7  for  Non  planar  Etching 

SAMPLE  release  1.7a  includes  programs  for  simulating  the  etching  of  nonplanar  layers  [Lyo88]. 
To  provide  an  interface  to  these  routines,  an  additional  command  was  added  to  SIMPL-2  to  create  a 
SAMPLE  input  file  for  nonplanar  etch  simulation.  To  begin,  the  top  of  the  profile,  where  the  topogra¬ 
phy  meets  air,  is  traversed  by  SIMPL-2  and  the  points  describing  this  contour  are  saved  in  an  array. 
The  polygon  containing  the  topmost  vertex  is  deleted,  unless  that  polygon  has  a  lowest  vertex  lower 
than  the  lowest  vertex  of  any  other  polygon  that  makes  up  the  surface.  Figure  3.4  shows  some  cases  of 
polygons  making  up  surface  layers,  and  which  would  be  deleted.  Once  a  polygon  is  deleted,  the  new 
surface  is  traversed  and  saved  as  the  next  layer.  Once  the  substrate  is  reached,  all  of  the  layers  are  writ¬ 
ten  to  a  SAMPLE  input  file  as  parameters  for  "nonplanar”  statements.  The  database  in  SIMPL-2  is 
restored  by  rebuilding  the  original  topography  from  data  which  had  been  stored  in  a  cross-section  data 
file. 

SIMPL-2  then  prompts  for  the  etch  rates  of  the  materials  represented  by  each  layer.  Finally, 
SIMPL-2  requests  the  etch  time  for  the  simulation.  Usually,  SIMPL  requests  the  etch  or  deposition  dis¬ 
tance  and  the  calculates  the  corresponding  time,  but  for  multiple  materials  with  different  etch  rates,  the 


resulting  etch  distance  is  known  only  after  the  simulation  is  run.  With  the  etch  rates  and  geometries  of 
each  of  the  layers,  the  interface  creates  a  complete  input  file  for  running  the  SAMPLE  etch  machine. 
Using  a  Unix  system  call,  SAMPLE  1.7a  is  run  with  the  input  file  and  the  results  are  stored  in  a  SAM¬ 
PLE  format  plot  file.  SIMPL-DIX  can  display  the  results  of  the  etch  simulation  in  the  profile  viewport 
as  shown  in  Figure  3.5.  Currently,  SIMPL-2  cannot  take  the  result  of  the  etch  simulation  and  update  its 
profile  data  base,  except  in  one  special  case:  if  the  etch  contour  consists  of  only  one  material  from  end 
to  end,  all  the  polygons  above  it  can  be  removed,  and  the  etch  contour  used  as  the  surface  of  a  new 
polygon  to  replace  the  polygon  currendy  cut  by  the  etch  contour.  This  is  certainly  not  a  general  inter¬ 
face,  but  it  can  be  useful  as  will  be  seen  in  a  later  example.  A  general  interface  will  be  available  in  the 
near  future,  but  is  not  a  part  of  this  project.  The  routines  used  for  the  nonpianar  etch  interface  are  listed 
in  Appendix  B  as  F77Layers.c.  These  routines  are  compiled  as  part  of  SIMPL-2. 


14 


4.  SIMPL  for  Analysis  of  a  Patterned  Photoresist  Planarization  Process 

4.1.  Process  Overview 

A  standard  technique  for  planarization  is  etch  back  of  a  resist-coated  dielectric  layer  [Ada81].  In 
this  technique  a  dielectric  film  is  deposited  over  the  topography  to  be  planarized.  An  organic  film  like 
photoresist  is  then  spun  on  in  such  a  way  as  to  planarize  the  entire  surface  of  the  layer.  The  combina¬ 
tion  organic  film  and  dielectric  are  then  etched  in  a  plasma  environment  that  has  been  designed  to  pro¬ 
duce  equal  etch  rates  in  both  materials.  The  limitations  of  this  process  have  been  documented  [Sti87], 
The  thickness  of  the  deposited  film  is  a  function  of  feature  density.  As  the  distance  between  the  features 
increases,  the  planarizing  effect  of  the  film  is  lost. 

A  planarization  process  using  a  sacrificial  fill  layer  of  patterned  photoresist  has  been  proposed 
which  solves  planarization  problems  encountered  in  both  intermetal  dielectric  for  a  1.2  micron  256K 
SRAM  technology  and  trench  isolation  for  a  0.8  micron  1M  SRAM  technology  [She88].  Photoresist  is 
used  to  fill  valleys  in  the  deposited  film.  After  the  dielectric  (e.g.  boro-silicate  glass)  is  deposited  on 
the  topography,  a  photoresist  is  spun  on  with  a  nominal  thickness  equal  to  the  step  height  of  the  under¬ 
lying  topography.  The  resist  is  the  patterned  in  such  a  way  as  to  remain  in  areas  where  conformal  cov¬ 
erage  of  the  dielectric  is  expected  to  occur.  A  second  layer  of  resist  is  now  deposited  on  top.  A  high 
degree  of  planarization  exists  for  the  last  deposition  step  since  most  of  the  nonplanar  regions  were  filled 
in  by  the  patterned  photoresist.  The  etch-back  proceeds  from  this  point  leaving  a  planar  dielectric  layer 
covering  the  topography. 

Misalignment  of  the  mask  used  to  pattern  the  first  layer  of  photoresist  can  be  expected  to  coun¬ 
teract  the  planarizing  benefits  of  this  process  and  degrade  device  performance.  However,  as  will  be 
shown  using  process  simulation,  this  masking  step  has  noncritical  dimensions  and  alignment  require¬ 
ments.  Adding  masking  steps  adds  to  the  cost  of  the  process,  but  some  of  this  cost  may  be  saved  given 
the  relatively  high  tolerance  of  this  masking  step  to  misalignment  error. 
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4.2.  SIMPL  Simulation  of  the  Patterned  Photoresist  Planarization  Process 

This  patterned  photoresist  planarization  process  can  be  simulated  using  SIMPL.  The  substrate  SI 
is  chosen  to  represent  some  underlying  topography.  In  this  case,  the  topography  consists  of  1.6/1. 4 
micron  lines  and  spaces  next  to  a  14  micron  region  with  no  lines.  The  height  of  the  lines  is  0.85 
microns.  Figure  4.1  shows  this  topography.  1.5  microns  of  glass  (PSG)  is  deposited  using  isotropic 
deposition.  0.75  microns  of  resist  (RST)  is  deposited  vertically  on  top  of  the  dielectric  as  shown  in  Fig¬ 
ure  4.2.  This  resist  is  patterned  with  the  mask  NB  leaving  a  box  of  resist  as  shown  in  Figure  4.3.  1.0 
microns  of  a  second  resist  is  deposited  as  shown  in  Figure  4.4.  After  using  SAMPLE  to  simulate  the 
etch  back,  the  resulting  profile  is  as  shown  in  Figure  4.5.  To  complete  the  process,  a  metal  layer  is 
deposited  on  top  of  the  dielectric  as  shown  in  Figure  4.6.  The  complete  SIMPL-2  process  file  for  this 
simulation  is  listed  in  Appendix  A  as  process.rbxl 

This  process  is  sensitive  to  misalignment  of  the  mask  used  to  pattern  the  resist  (NB).  The 
WORST  feature  of  SIMPL-DIX  is  used  to  shift  the  mask  NB  to  the  right.  The  result  of  shifting  the 
mask  NB  to  the  right  by  0.6,  0.9,  1.2,  1.8  microns  and  then  simulating  the  process  is  shown  in  Figures 
4.7,  4.8,  4.9,  and  4.10  respectively.  Also,  if  the  mask  NB  is  originally  too  small  due  to  excessive  pro¬ 
cess  bias,  the  result  is  shown  in  Figure  4.11.  At  a  certain  point,  the  effect  of  mask  misalignment  is  to 
greatly  diminish  the  planarizing  advantage  of  the  process.  After  0.7  microns  misalignment,  the  prob¬ 
lems  are  severe. 

4.3.  Results  -  RACPLE  Analysis 

To  get  a  more  quantitative  measurement  of  the  effect  of  mask  misalignment  in  this  process,  RAC¬ 
PLE  can  be  used  to  measure  the  number  of  lateral  squares  in  the  dielectric  layer.  Figure  4.12  is  a  plot 
of  lateral  squares  of  PSG  versus  misalignment.  For  no  misalignment  and  misalignment  of  0.6  microns, 
RACPLE  measures  about  8.2  lateral  squares  for  the  PSG  layer.  At  0.7  microns  misalignment,  the  RAC¬ 
PLE  measurement  jumps  to  12.4  and  continues  to  climb  to  13.9  for  a  misalignment  of  0.9  microns.  For 
1.1  microns  misalignment  RACPLE  measures  12.28  squares,  and  for  more  than  1.2  microns  misalign¬ 
ment,  the  measurement  settles  at  11.33  lateral  squares  and  remains  there  for  misalignment  of  up  to  2.4 
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microns.  Similarly,  for  the  profile  of  Figure  4.11,  in  which  the  mask  NB  is  biased  too  small,  RACPLE 
calculates  1 1 .83  lateral  squares  for  the  dielectric  layer. 

It  is  apparent  that  the  misalignment  or  mask  size  error  to  create  a  noticeable  degradation  of  dev¬ 
ice  performance  is  on  the  order  of  50%  of  the  minimum  linewidth  for  this  example.  This  supports  the 
claim  that  the  masking  step  has  noncritical  dimensions.  The  degree  to  which  the  masking  step  tolerates 
error  has  been  shown  here  with  SIMPL.  One  interesting  feature  of  Figure  4.12  is  the  sudden  jump  in 
the  number  of  lateral  squares  at  0.7  microns  misalignment  which  settles  at  around  1.2  microns  misalign¬ 
ment.  One  cause  of  this  effect  is  the  fact  that  at  above  1.2  microns  misalignment,  the  thinning  of  the 
dielectric  on  the  left  side  of  the  structure  is  compensated  by  an  excess  of  dielectric  on  the  right  side.  At 
0.7  to  0.9  microns  misalignment,  there  is  a  serious  thinning  on  the  left  side,  but  no  noticeable  thicken¬ 
ing  on  the  right.  The  capacitance  across  the  dielectric  on  the  left  side  is  actually  quite  high  for  all 
cases  of  misalignment  greater  than  50%  of  the  minimum  linewidth,  but  is  compensated  on  the  right  side 
if  the  mask  is  misaligned  enough. 

Using  SIMPL,  it  has  been  shown  that  mask  misalignment  tolerance  is  good  up  to  about  0.6 
microns  for  this  case,  but  results  in  serious  performance  degradation  for  further  misalignment.  The 
effects  of  altering  layer  thicknesses,  or  changing  feature  dimensions  and  spacings  can  be  expected  to 
affect  the  misalignment  tolerance.  These  effects  can  be  investigated  with  SIMPL. 

4.4.  Comments  on  Approach 

This  approach  to  analyzing  the  above  planarization  process  can  identify  certain  trends  using  ele¬ 
mentary  models.  It  would  be  interesting  to  include  additional  process  effects  and  more  comprehensive 
analysis  of  the  exact  electrical  nature  of  the  device.  Microloading  effects  which  would  cause  a  variation 
in  etch  rate  across  an  individual  die  are  not  accounted  for  by  the  etch  simulation  models.  As  a  result, 
certain  topography  effects  which  might  be  problematic  in  real  devices  are  not  seen  in  this  simulation. 
Additionally,  RACPLE  gives  only  an  estimate  of  the  electrical  properties  of  the  profile.  A  complete 
analysis  of  the  multiple  parasitic  capacitances  and  even  inductances  would  be  of  interest. 


17 


5.  Application  to  VLSI  Hopfield  Neural  Networks 

5.1.  SIMPL  for  Neural  Network  Analysis 

It  was  originally  proposed  that  the  linking  of  SIMPL  process  simulation  with  RACPLE  analysis 
could  be  used  to  investigate  parasitic  resistive  and  capacitative  loading  effects  in  highly  interconnected 
computational  structures  such  as  neural  networks.  SIMPL  and  RACPLE  have  been  applied  to  an 
analysis  of  parasitic  loading  in  a  structure  similar  to  that  developed  at  AT&T  Bell  Laboratories  [Jac86]. 
It  will  be  shown  that  layout  and  topographical  features  in  processing  do  have  an  impact  on  neural  net¬ 
work  performance. 

5.2.  Overview  of  a  Reported  VLSI  Neural  Network 

The  properties  of  highly  interconnected  arrays  of  amplifiers  have  generated  much  interest  for  their 
potential  use  in  a  new  class  of  computing  circuits.  The  properties  of  such  networks  have  much  in  com¬ 
mon  with  biological  information  processing  systems  (brains)  in  that  they  are  massively  parallel  and  fault 
tolerant.  The  basic  network  configuration  is  shown  at  the  top  of  Figure  5.1.  Several  amplifiers  are  con¬ 
nected  such  that  each  amplifier  output  is  available  as  input  to  any  of  the  other  amplifiers.  The  actual 
feedback  connection  is  made  with  a  resistor,  and  the  pattern  of  resistors  in  the  network  determines  the 
behavior  of  the  entire  circuit. 

The  basic  operation  of  the  circuit  can  be  described  as  analogous  to  the  motion  of  a  particle 
through  a  potential  energy  field  in  a  multidimensional  space.  This  analogy  holds  if  the  matrix  of  inter¬ 
connect  resistances  is  symmetric.  The  output  voltages  V ,  of  each  of  the  N  amplifiers  are  independent 
coordinates  in  space.  The  amplifier  gain  characteristics  are  assumed  to  be  symmetric  around  V=0.  For 
this  case,  the  energy  function  is 

E=I?tjViVj+±\g-'<y)dv 
<</  r‘ 

where  V=g(u)  is  the  transfer  function  of  the  amplifiers,  T,j=T }l=\!RL]  are  the  coupling  resistors 
between  the  amplifiers,  1/-/?,=^  1//?,,  and  T„= 0  [Hop84],  If  the  system  is  put  into  any  particular  state 
by  applying  voltages  at  the  inputs,  the  energy  function  gradient  w'ill  cause  the  circuit  to  relax  to  a  stable 
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state  which  is  close  to  the  initial  state.  This  type  of  circuit  can  be  used  as  a  content-addressable 
memory  [How87],  The  search  operation  is  done  in  a  fully  parallel  way  and  the  time  to  reach  a  solution 
is  determined  by  the  speed  of  the  amplifiers  and  the  time  constant  of  the  resistor  network. 

At  the  bottom  of  Figure  5.1  is  a  circuit  diagram  showing  the  basic  make-up  of  a  single  neuron. 
The  resistances  and  capacitances,  in  addition  to  the  resistive  weight,  are  due  to  parasitics  in  the  VLSI 
implementation  of  the  network.  To  study  these  effects,  a  single  interconnect  element  was  fabricated 
with  SIMPL. 


5-3.  SIMPL  Simulation  of  a  Neural  Network  Element 

Using  electron-beam  lithography,  a  12x12  resistor  matrix  that  fits  into  a  6x6  micron  square  was 
fabricated  by  a  group  at  AT&T  Bell  Laboratories  [Jac86].  To  generate  a  device  cross  section  using 
SIMPL,  the  basic  geometry  of  the  AT&T  device  was  used.  A  0.5  micron  thick  oxide  layer  was  grown 
on  the  substrate.  0.1  microns  of  tungsten  was  deposited  and  patterned  into  0.3  micron  lines  and  spaces. 
0.1  microns  of  a  polyimide  was  deposited  as  an  interplanar  dielectric  (using  SIMPL  anisotropic  deposi¬ 
tion  with  a  20  degree  source  angle.  A  resistor  hole  was  created  and  filled  with  polysilicon.  Finally,  a 
0.1  micron  layer  of  nickel  was  deposited  (using  SIMPL  vertical  deposition).  The  resulting  cross  section 
is  shown  in  Figure  5.2.  The  cut  line  was  positioned  such  that  the  profile  generated  represents  one  period 
(0.5  microns)  of  the  tungsten  line.  A  second  simulation  for  a  cross-over  point  with  no  contact  is  shown 
in  Figure  5.3.  The  SIMPL  process  file  for  Figure  5.3  is  listed  in  Appendix  A. 


5.4.  Parasitic  Analysis 

For  a  planar  nickel  layer  which  is  0.5  microns  long,  0.1  microns  thick  and  0.2  microns  wide.  The 
number  of  lateral  squares  is  5  and  the  overall  resistance  is 


R  =7.85£I-cmx5x 

The  inter-metal  capacitance  is  on  the  order  of 


1 

0.2|im 


C  =3.9x8. 85x  10'14— x3x0.2pjn 
cm 

Using  RACPLE  to  measure  the  resistance  of  the  nickel  electrode,  the  number  of  lateral  squares  is  7.085 
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instead  of  5,  giving  an  increase  of  40%.  Measuring  the  inter  metal  capacitance  by  running  RACPLE 
over  the  width  of  the  tungsten  electrode  only,  gives  4.7  lateral  squares  instead  of  3,  an  increase  of  more 
than  56%.  Additional  capacitances  from  the  electrode  to  the  substrate  are  negligible  in  comparison  with 
the  inter-metal  capacitance.  Using  these  values  in  the  above  equations,  the  total  RC  of  the  upper  elec¬ 
trode  is  about  1.0  picoseconds.  If  there  are  a  thousand  such  tungsten  lines,  the  total  RC  delay  of  the 
nickel  line  is  about  1  microsecond.  If  we  assume  an  amplifier  delay  on  the  order  of  a  microsecond, 
clearly  the  interconnect  delay  time  is  an  important  mechanism.  Topography  related  effects  can  have  a 
serious  impact  on  neural  network  circuit  performance,  since  the  total  distributed  RC  delay  constant  of 
the  counter  electrode  increases  as  the  square  of  the  number  of  neurons.  With  planarization  techniques, 
the  interconnect  delay  can  be  reduced  to  improve  circuit  performance. 

5.5.  Comments  on  this  Approach 

SIMPL  and  RACPLE  are  in  place  for  studying  topography  and  process  dependent  effects  in  novel 
circuits  such  as  neural  networks.  As  demonstrated  above,  basic  trends  in  the  characteristics  of  novel 
devices  can  be  investigated  with  integrated  process  and  device  simulation.  It  is  worth  pointing  out  that 
industry  is  not  yet  developing  VLSI  neural  networks  which  use  processes  much  different  from  those 
used  for  conventional  CMOS  chips.  The  interest  in  neural  networks  now  is  in  getting  the  chips  to  work 
and  finding  applications  for  them  [Hec88].  Still,  many  new  architectures  and  designs  for  VLSI  neural 
networks  are  being  proposed,  and  tools  such  as  SIMPL  and  RACPLE  can  be  used  to  provide  an  initial 
assessment  of  some  of  the  electrical  issues  involved. 

It  would  be  interesting  to  perform  a  full  three-dimensional  analysis  of  the  interconnect  structure, 
including  ones  with  resistive  weights.  Process  simulation  should  generate  the  device  geometry  in  three 
dimensions  to  provide  the  right  link  between  process  flow  and  device  analysis.  The  nature  of  the  circuit 
poses  some  interesting  challenges  as  well.  It  has  been  proposed  that  the  stability  properties  of  Hopfield 
neural  networks  can  be  related  to  the  properties  of  individual  neurons  [Mic87].  This  synthesis  of  system 
theory  and  device  technology  should  also  be  brought  to  the  analysis.  Rigorous  software  for  process  and 
device  simulation  will  be  needed  for  problems  of  dense  interconnect  networks  in  general,  and  has  a 


ready  appiication  in  neural  networks. 
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6.  Future  Directions  for  SIMPL 

6.1.  Profile  Data  Management 

In  the  area  of  profile  data  management,  two  questions  must  be  addressed.  First,  what  data  should 
be  maintained?  Second,  how  should  this  data  be  transferred  to  places  where  it  is  needed?  For  SIMPL, 
the  first  question  can  be  answered  this  way:  SIMPL  should  maintain  a  description  of  a  device  profile 
which  completely  describes  the  device  geometry,  and  materials  involved.  The  profile  data  description  as 
it  now  exists  handles  most  of  this  information.  However,  some  additions  are  needed  to  completely 
describe  a  profile.  Currently,  SIMPL  only  describes  the  net  charge  of  dopant  contained  at  a  particular 
location  in  the  substrate.  This  should  be  changed  so  that  the  species  of  dopant  is  included.  Such  infor¬ 
mation  is  needed  for  diffusion  simulation  with  programs  like  SUPREM.  SIMPL  should  also  maintain 
information  about  impurities  in  materials  other  than  the  silicon  substrate.  No  information  is  currently 
maintained  about  impurities  present  in  a  gate  oxide,  for  example.  Another  addition  to  the  profile 
description  which  is  still  needed  is  the  ability  to  describe  a  floating  island  of  material  which  is  com¬ 
pletely  surrounded  by  another  material.  There  are  ways  around  this,  usually  involving  dividing  up  the 
surrounding  material  into  two  parts.  Ideally  however,  the  SIMPL  data  base  would  handle  a  floating 
island  of  material  as  a  normal  case. 

The  second  question  involves  how  data  is  transferred.  Currently,  this  is  one  of  the  major 
bottlenecks  which  noticeably  slows  down  the  performance  of  SIMPL.  SIMPL-DIX  and  SIMPL-2  use 
files  to  transfer  profile  information  back  and  forth.  The  time  spent  writing  files  can  be  significant  when 
large  pieces  of  data  must  be  transferred.  It  is  easy  to  write  programs  that  communicate  through  files,  but 
it  is  not  very  efficient  from  a  user’s  point  of  view.  In  this  day  of  powerful  computer  work-stations,  with 
several  megabytes  of  RAM  storage  available,  SIMPL  should  avoid  communication  with  files.  Instead, 
data  should  be  stored  in  memory  common  to  both  SIMPL-DIX  and  SIMPL-2. 

Likewise,  in  the  area  of  communication  with  other  programs,  it  is  desired  to  get  around  the  need 
to  use  files.  This  is  not  an  easy  problem  to  solve  yet,  especially  for  programs  written  in  different  pro¬ 
gramming  languages.  There  is  some  promise,  however,  that  a  PIF  data  base  will  one  day  become 
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available.  Eventually  it  will  be  possible  to  send  a  program  a  pointer  to  the  root  of  a  data  tree  instead  of 
a  file  full  of  data.  Without  a  PIF  parser  this  is  not  yet  a  possible  alternative,  but  it  should  be  pursued. 

6.2.  Tools  Integration  with  SIMPL 

There  is  already  strong  interest  in  industry  for  an  integrated  design  environment  for  process  and 
device  simulation.  SIMPL  has  the  potential  for  filling  this  gap,  but  work  remains  to  be  done.  SIMPL- 
DIX  is  a  suitable  framework  for  creating  menu  driven  interfaces  to  various  programs,  but  it  is  not 
always  an  easy  environment  in  which  to  do  development  work.  For  a  programmer  wishing  to  integrate  a 
new  simulator  into  SIMPL-DIX,  it  would  be  useful  to  have  some  sort  of  library  of  standardized  routines 
from  which  to  build  the  application  interface.  Many  routines  which  fit  that  description  already  exist  in 
SIMPL -2  and  SIMPL-DIX  but  they  are  not  organized  in  an  efficient  manner.  Many  must  also  be  rewrit¬ 
ten  in  order  to  be  useful  for  general  applications.  As  it  stands  now,  it  often  takes  programmers  several 
months  to  understand  the  intricacies  of  SIMPL.  Instead  of  tracking  down  routines  which  already  exist,  it 
is  often  easier  to  write  new  ones.  This  only  adds  to  the  size  and  complexity  of  SIMPL.  One  might  say 
that  SIMPL  is  becoming  so  convoluted  that  the  name  should  be  changed  to  HARD.  SIMPL,  if  fully 
integrated  with  programs  like  SAMPLE,  SUPREM  and  PISCES,  would  become  a  powerful  design 
environment.  With  a  standardized  approach  to  integrating  new  simulation  software  into  its  framework, 
SIMPL  would  remain  at  the  leading  edge  in  the  face  of  rapidly  changing  technology.  Admittedly,  this 
would  be  a  major  undertaking,  but  industry  would  be  very  interested  in  a  design  environment  which 
could  easily  be  connected  to  proprietary  device  and  process  simulation  software. 

Listed  in  Appendix  C  are  several  of  the  C  functions  in  SIMPL-2  which  are  useful  for  tools 
integration  and  profile  manipulation.  This  list  includes  C  functions  available  in  the  most  recent  release 
version  of  SIMPL-2  as  well  as  new  ones  written  as  a  part  of  the  ongoing  development  of  SIMPL.  A 
few  routines  in  SIMPL-DIX  which  were  used  in  this  project  are  also  listed.  The  routines  are  organized 
according  to  their  use  in  tools  integration  and  profile  manipulation.  It  is  not  an  exhaustive  list  of  the 
functions  available  in  SIMPL  but  is  intended  as  a  first  step  in  organizing  the  routines  in  a  useful 


fashion. 
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Another  area  to  be  considered  is  the  possibility  of  integrating  SIMPL  with  the  OCT/VEM 
environment  [Har86].  Many  of  the  software  tools  needed  for  integrating  SIMPL  with  other  simulators 
exists  in  OCT/VEM.  Additionally,  the  issue  of  upgrading  from  X  windows  version  10  to  version  1 1 
would  automatically  be  solved  if  SIMPL  is  included  in  OCT/VEM.  Once  those  programs  switch  to  XI 1, 
so  will  SIMPL.  The  trade-offs  and  advantages  of  integrating  SIMPL  with  OCT/VEM  should  be  further 
explored.  Some  combination  of  these  options,  where  SIMPL  can  be  used  as  a  stand-alone  environment 
and  where  it  is  also  available  in  OCT/VEM  may  be  the  most  flexible  and  useful  approach. 

A  standard  PEF  format  is  also  a  necessity  in  order  to  pursue  the  goal  of  tools  integration.  It  is  pos¬ 
sible  to  develop  custom  interfaces  for  each  simulation  program  that  becomes  available,  but  the  time 
involved  is  often  great  If  every  process  and  device  simulation  program  communicated  with  PIF,  this 
aspect  of  tools  integration  would  be  trivial. 

6 3.  Problems  in  the  Implementation  of  SIMPL 

Recent  releases  of  SIMPL  suffer  from  some  implementation  problems  which  hamper  the 
effectiveness  of  the  program.  The  biggest  problems  involve  the  internal  etch  simulation  models.  Many 
standard  cases  are  not  handled  correctly,  and  the  resulting  profiles  are  incorrect.  There  are  also  cases 
where  SIMPL-2  incorrectly  searches  the  database,  which  instead  of  merely  giving  a  false  profile,  cause 
the  program  to  crash  unexpectedly.  The  implantation  routines  often  suffer  from  overflow  problems. 
The  grid  which  SIMPL-2  uses  to  store  doping  concentration  information  is  sometimes  allocated 
incorrectly.  The  routines  are  also  prone  towards  developing  infinite  loops  which  freeze  the  program. 
SIMPL-DIX  has  some  minor  problems,  generally  resulting  when  some  global  variable  is  not  reset  prop¬ 
erly.  Occasionally,  attempts  to  re-initialize  SIMPL-2  from  SIMPL-DIX  are  not  successful. 

This  lack  of  robustness  takes  away  from  the  usefulness  and  credibility  of  the  program.  When  a 
perfectly  normal  process  cannot  be  simulated  by  SIMPL,  most  potential  users  simply  give  up.  These 
problems  are  not  due  to  major  shortcomings  of  the  ideas  behind  SIMPL,  but  can  be  traced  to  problems 
of  implementation  and  inconsistencies  in  the  software.  At  some  point,  these  implementation  problems 
will  need  to  be  thoroughly  investigated  and  cleaned  up. 
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6.4.  Conclusion 

SIMPL  is  a  useful  tool  for  studying  the  complex  interrelationship  of  physical  phenomena  that  go 
into  modem  integrated  circuit  process  design.  This  project  has  demonstrated  how  some  basic  tools  can 
be  integrated  into  the  framework  of  SIMPL,  shown  how  SIMPL  can  be  used  to  analyze  proposed 
processes,  and  listed  some  of  the  work  that  must  be  done  to  keep  SIMPL  at  the  cutting  edge  of  CAD 
technology  for  process  and  device  simulation.  It  is  hoped  that  this  report  is  a  useful  contribution  to  this 
increasingly  important  field  in  integrated  circuit  design  and  manufacture. 
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Figures 


Figure  1.1.  SIMPL  as  an  all-purpose  design  interface 


vertex  d 


i  =  vertex 


Figure  2. 1 .  Linked  polygon  data  structure 


Figure  2.2.  Polygons  represented  as  layers 


Figure  2.3  Example  of  displaying  RACPLE  data  with  SEMPL 


Figure  3.1.  typical  polygon  with  vertical  edges. 


1 


Figure  3.2.  head  vertex  classifications. 


Figure  3.3.  polygon  split  into  two  layers. 


Case  1 .  One  polygon  with  top 
vertex  equal  to  highest  point 
on  surface. 


Case2.  Polygon  nested  in  a  valley 
created  by  another  polygon. 


Case  3.  Multiple  polygons  with 
vertexes  equal  to  highest  point 
in  surface. 


Polygon  identified  as  top  polygon. 


Figure  3.4.  Examples  identifying  which  surface  polygon 
would  be  deleted  in  the  routines  to  write  SIMPL-2  profiles 
as  a  series  of  layers. 


Figure  3.5  Example  of  displaying  SAMPLE  non-planar  etch  with  SIMPL 


Figure  4.1  Underlying  topography,  1.6/1. 4  micron  lines  and  spaces 


Figure  4.2  After  depositing  dielectric  and  first  photoresist. . . 


Figure  4.3  After  patterning  photoresist.  . . 


Figure  4.4  After  depositing  second  layer  of  photoresist. . . 


Figure  4.5  Topography  remaining  after  etch-back 


Figure  4.6  Metallization  step  gives  final  topography. 


Figure  4.7  Effect  of  0.6  micron  misalignment 


Figure  4.8  Effect  of  0.9  micron  misalignment 


Figure  4.9  Effect  of  1.2  micron  misalignment 


Figure  4. 10  Effect  of  1 .8  micron  misalignment 


Figure  4.1 1  Effect  of  mask  NB  being  too  small 


Figure  4.12.  Effective  lateral  squares  in  the  dielectric 
versus  misalignment  of  mask  used  to  pattern 
the  resist  in  the  patterned  resist  planarization 
process. 


Figure  5.2.  One  interconnect  element  of  a  VLSI  Hopfield  Neural  network. 

In  this  example,  POLY  is  polyimide,  PSG  is  polysilicon,  WMTL 

is  tungsten  and  NMTL  is  nickel.  The  thickness  of  each  layer 
is  approximately  0.1  microns.  The  width  of  the  tungsten  lines 
is  0.3  microns,  and  their  spacing  is  0.2  microns. 


Figure  5.3.  Counter  electrode  crossing  base  electrode  in  VLSI  Hopfield 
neural  network  with  no  resistive  contact. 


Appendix  A 


SIMPL-2  Process  Files 


*  *  ★  ***  ★  * *★★*★***★★*****★*★*★****★*★**★★★★★★★★★*★*★*★★**★★★****★★★★★★★★★*★*★★★★★'** 
LAYOUT  FILE  :  rbxl.cif 

SUBSTRATE  TYPE:  p  type,  concentration  le!4 

CUT-LINE  COORDINATES  :  xl  -  -1705,  yl  -  271 

x2  =  596,  y2  ■=  276 

★  ★★★★★★****************************************************************1'"*"'***''’"' 
*  i  * 

WHICH  PROCESS  ?  DEPO 
NAME  OF  THE  MATERIAL  ?  rst 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  1.0 

ISO,  ANISO,  VERT,  or  SAMPLE  MENU  (I,A,V,or  S)  ?  v 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  2  * 

WHICH  PROCESS  ?  EXPO 
WHICH  MASK  ?  ns 

INVERT  THE  MASK  (yes  or  no)  ?  no 
NAME  OF  THE  EXPOSED  RESIST  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  3  * 

WHICH  PROCESS  ?  DEVL 

NAME  OF  THE  LAYER  TO  BE  DEVELOPED  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  4  ★ 

WHICH  PROCESS  ?  ETCH 

WHICH  LAYER  DO  YOU  WANT  ETCH  ?  si 

ETCH  ALL  (yes  or  no)  ?  no 

AMOUNT  OF  VERTICAL  ETCH  (micro_meter )  ?  0.85 

RATIO  X/Z  OF  ETCHING  (0.0  <=  RATIO  <=  1 . 0 )  ?  0.0 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  5  * 

WHICH  PROCESS  ?  DEVL 

NAME  OF  THE  LAYER  TO  BE  DEVELOPED  ?  rst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  6  * 

WHICH  PROCESS  ?  DEPO 
NAME  OF  THE  MATERIAL  ?  psg 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  1.5 

ISO,  ANISO,  VERT,  or  SAMPLE  MENU  (I,A,V,or  S)  ?  i 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


★  7  ★ 


WHICH  PROCESS  ?  DEPO 


NAME  OF  THE  MATERIAL  ?  rst 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  0.75 

ISO,  ANISO,  VERT,  or  SAMPLE  MENU  (I,A,V,or  S)  ?  v 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  8  * 

WHICH  PROCESS  ?  EXPO 
WHICH  MASK  ?  nb 

INVERT  -THE  MASK  (yes  or  no)  ?  no 
NAME  OF  THE  EXPOSED  RESIST  ?  erst 

DO  YOU -WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  9  * 

WHICH  PROCESS  ?  DEVL 

NAME  OF  THE  LAYER  TO  BE  DEVELOPED  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  10  * 

WHICH  PROCESS  ?  DEPO 
NAME  OF  THE  MATERIAL  ?  rst2 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  1.0 

ISO,  ANISO,  VERT,  or  SAMPLE  MENU  (I,A,V,or  S)  ?  i 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


★  11  * 

WHICH  PROCESS  ?  ETCN 

etchrate  for  RST2,  layer  3  (um/sec)  ?  0.01 
etchrate  for  RST,  layer  2  (um/sec)  ?  0.01 
etchrate  for  PSG,  layer  1  (um/sec)  ?  0.01 
etchrate  for  SI,  layer  0  (um/sec)  ?  0.001 
timestep  in  seconds  ?  50 
number  of  steps  ?  4 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  12  * 

WHICH  PROCESS  ?  ETCU 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


*  13  * 

WHICH  PROCESS  ?  DEPO 

NAME  OF  THE  MATERIAL  ?  metl 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  0.85 

ISO,  ANISO,  VERT,  or  SAMPLE  MENU  (I,A,V,or  S)  ?  i 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


WHICH  PROCESS  ?  END 


LAYOUT  FILE  :  nn.cif 

SUBSTRATE  TYPE:  p  type,  concentration  0 

CUT-LINE  COORDINATES  :  xl  -  -59,  yl  -  12 

x2  -  30,  y2  -  12 


*  1  * 

WHICH  PROCESS  ?  OXID 

OXIDE  THICKNESS  (micro-meter)  ?  .5 

Xt  (micro-meter)  ?  .5 

Xe  (micro-meter)  ?  .25 

ul  ?  .1 

u2  ?  .5 

u3  ?  .9 

dl  ?  .1 

d2  ?  .5 

d3  ?  .  9 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 
WHICH  PROCESS  ?  WAIT 
*  2  * 

WHICH  PROCESS  ?  DEPO 
NAME  OF  THE  MATERIAL  ?  wmtl 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  .1 
ISOTROPIC,  ANISOTROPIC,  OR  VERTICAL  (I,  A,  or  V)  ?  v 
DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 

WHICH  PROCESS  ?  WAIT 

*  3  * 

WHICH  PROCESS  ?  DEPO 
NAME  OF  THE  MATERIAL  ?  rst 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  .25 
ISOTROPIC,  ANISOTROPIC,  OR  VERTICAL  (I,  A,  or  V)  ?  v 
DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 

WHICH  PROCESS  ?  WAIT 

*  ^  ★ 

WHICH  PROCESS  ?  EXPO 

WHICH  MASK  ?  wmtl 

INVERT  THE  MASK  (yes  or  no)  ?  no 

NAME  OF  THE  EXPOSED  RESIST  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 
WHICH  PROCESS  ?  WAIT 
*  5  * 

WHICH  PROCESS  ?  DEVL 

NAME  OF  THE  LAYER  TO  BE  DEVELOPED  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ?  yes 


WHICH  PROCESS  ?  WAIT 


*  6  * 


WHICH  PROCESS  ?  ETCH 

WHICH  LAYER  DO  YOU  WANT  ETCH  ?  wmtl 
ETCH  ALL  (yes  or  no)  ?  no 

AMOUNT  OF  VERTICAL  ETCH  (micro_meter )  ?  .1 
RATIO  X/Z  OF  ETCHING  (0.0  <-  RATIO  <-  1 . 0 )  ?  .1 
DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ' 

WHICH  PROCESS  ?  WAIT 

*  7  * 

WHICH  PROCESS  ?  DEVL 

NAME  OF  THE  LAYER  TO  BE  DEVELOPED  ?  rst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no)  ' 

WHICH  PROCESS  ?  WAIT 

*  8  * 

WHICH  PROCESS  ?  DEPO 

NAME  OF  THE  MATERIAL  ?  poly 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  .10 
ISOTROPIC,  ANISOTROPIC,  OR  VERTICAL  (I,  A,  or  V)  ? 
SPUTTERING  SOURCE  ANGLE  (degrees)  ?  20 
DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no) 

WHICH  PROCESS  ?  WAIT 

*  12  * 

WHICH  PROCESS  ?  EXPO 
WHICH  MASK  ?  psg 

INVERT  THE  MASK  (yes  or  no)  ?  yes 
NAME  OF  THE  EXPOSED  RESIST  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no) 
WHICH  PROCESS  ?  WAIT 

*  13  * 

WHICH  PROCESS  ?  DEVL 

NAME  OF  THE  LAYER  TO  BE  DEVELOPED  ?  erst 

DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no) 

WHICH  PROCESS  ?  WAIT 

*  14  * 

WHICH  PROCESS  ?  DEPO 
NAME  OF  THE  MATERIAL  ?  nmtl 

THICKNESS  OF  THE  MATERIAL  (micro-meter)  ?  .10 
ISOTROPIC,  ANISOTROPIC,  OR  VERTICAL  (I,  A,  or  V)  ? 
DO  YOU  WANT  TO  DRAW  THE  CROSS  SECTION  (yes  or  no) 

WHICH  PROCESS  ?  WAIT 


yes 


:  yes 


a 

?  yes 


?  yes 


?  yes 


v 

?  yes 


WHICH  PROCESS  ?  END 


Appendix  B 


Source  code  for  the  majority  of  the  routines  developed  for  this  project. 


dix_actions4.c 

F77Layers.c 


dix  action4 . c 


Thu  Dec  1  15:54:34  1988 


1 


/*  dl *  act  1 on4 . c 

*  Fourth  part  of  SIHPL-DIX  action  routines. 

• 

•  See  "comma ndcont rol .c" . 

* 

•  Edward  W.  Scheckler  Nov.  15,  1988 

*  Copyright  (C)  1908  U.  C.  Berkeley  SAMPLE  Group 

*/ 

•Include  <stdlo.h> 

•Include  <math.h> 

•Include  <X/Xlib.h> 

•Include  "slmpl-dix.h" 

•Include  "slmpl.h" 

•Include  "cif.h" 

•Include  "command. h" 

•Include  "display. h" 

•Include  "default. h" 


/* 

*  External  declarations. 

*/ 

char  *RACPLE  Path; 

extern  dlxVlewport  DIX  Viewport  (VIEWPORT  SUE |  ; 

extern  floatVlew  Prof  1 le  Vlew; 
extern  floatBound  Profile  Bound; 


extern  short  PatternSl ze; 

extern  short  Commandld; 
extern  short  Menujd; 

extern  short  Cut  1 1 ne  Stat us; 
extern  short  Layout  Status; 
extern  short  PatternStat us; 
extern  short  Prof  1 le_St at us; 
extern  short  Prof  1 leFrame  Status; 


extern  Fontlnfo  • BodyFont Info; 
extern  Color  BackgroundColor; 
extern  Color  ForegroundColor; 
extern  slmplPolygon  •SIMPLPolygonRt ; 
static  float  x 1 ef t ( 50 | , y 1 ef t [ 50 | ; 
static  float  squares(50|; 

/*'***'* 

Routines  to  call  a  version  of  RACP1.E 
Added:  Oct  ?4,  1908  EWS 

If  the  file  RACPLF.  RESULTS  already  exists,  then  DIX 
assumes  that  It  belongs  to  the  current  profile.  If  not 
such  a  file  Is  created  by  calling  RACPLF.  for  each  polygon 
I  n  t  he  prof  1 1 e. 


/ 


DoRacple ( ) 

I 

FILE  *fpp,  *fopen{); 
slmplPolygon  ‘Polygonptr; 
char  polyname (NAME  SIZE  1; 
char  junk(80| , colon; 
struct  float  path  *Polypath; 

Int  1,11; 

Polygonptr  •  SIMPLPo lygonRt ; 

1  -  0; 

1 f ( (fpp-fopen ("RACPLE ^RESULTS", "r") 1 ! -NULL) 1 
1  f  (GetYesOrNo  ("  RACPLERESULTS  exists.  Use  It?")—  YES)  1 
/•read  the  file*/ 

f scanf < fpp, "Is  Is", junk, junk) ; 

f scanf ( fpp, "Is  Is  Is  Is  Is  Is  Is  Is  Is" , junk, junk, 
junk,  junk, junk,  Junk,  Junk,  junk, Junk) ; 

for  (;;)  I 

f scanf 1 fpp, "Id", ill ) ; 

If (11  ! -  1  >  break; 
f scanf ( fpp, "Is", polyname) ; 
f scanf (fpp, "If", ixleft l 1 1 ) ; 
f scanf ( fpp, "lf*,4yleft|l)); 

f scanf (fpp, “If", 4  squares! 1 )) ; 

1  ; 

I 

fclose ( fpp) ; 

)  else  \ 
fclose ( fpp) ; 

goto  makenew; 

1  /‘Get  YesOrNo*/ 

)  else  ( 

makenew:  fpp  -  fopen("RACPLERESULTS","w"); 

fprlntf (fpp, "•*********•**••• *RACPLE  RESULTS* *********** . \n") . 

fprlntf (fpp, "Material  Type,  Head  Vertex  (x, z) ,  •  of  Lateral  SquaresXn") ; 
while  (Polygonptr  f-  NULL)  I 

CreateRacplel nput (l , Polygonpt  r, ixleft { 1 ) ,  4y left ( 1 ) ) ; 

RunRacple (1 , 4  squares! 11); 
sprlntf (polyname, "Is", Pol ygonpt r->name) ; 

fprlntf (fpp, "13d  16s  112. 3f  17. 3f  I9.3f\n*,l, 

polyname, xleft ( 1 ) ,yleft ( 1 ),  squares! 1 1 ) ; 

Polygonptr  -  Po lygonpt r->next ; 
l  +  t; 

)  /‘while*/ 
fclose ( fpp) ; 

)  /"if*/ 

Print  RacpleData ( l ) ; 

I 

/******* 

Prl nt RacpleData (number  of  points) 

Routine  to  sort  each  row  of  dat a  by  value  of 
x  coordinate  and  print  It  to  profile  viewport. 

Revised  11/8/88  EWS 

*•*•***/ 

P  r  I  nt  Racp  leOat  a  ( 1  ) 
l  nt  I  ; 

( 

d  1  xVl  ewf>ort  view; 

Int  max  per  row, no  rows, no  col; 

Int  left  edge, right  edge,  top  edge, bottom  edge; 


dix  action4 . c 
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lnt  wval  ,  hval,  xval  (SO) ,yval  (50),  xpolnt  (50)  , ypolnt  I  50  J  ; 

int  Index, j, k,  1, lm, kin, gap,n, flag; 
float  tempx,  tempy,  temps, prevx; 

int  samexl  ( 1 0  J ,  samexr  ( 10 ) ,  nosane,  ltoh(10| , htol 110); 
char  value (151; 

char  promptstrlng(SO); 

Pixel  fore, back; 

Fontlnfo  ‘font; 

FILE  *fp; 

view  -  DIX  Viewport ( 1) ; 

left_edge  -  Get Prof 1 leVlewX (view,Profl leView, Prof  1 leBound. left) ; 
rlghtedge  -  GetProfl  leVlewX (view, Prof lie  View, Prof  1 le  Bound. right) ; 
topedge  -  GetProfl leVlewY (view, Prof  1 leView, Prof  1 le  Bound. top) ; 

bottomedge-GetProf  HeVlewY  (view.  Prof  1  leView, Prof 1 leBound. bottom)  ; 
sprint  f (value, "%f ",  squares (0J ) ; 
wval  -  GetTextWldth  (value,  5,  Body  Font  Info)  ; 
hval  -  (2  *TF.XT_MARGIN) ♦Body  Font  I nfo->height; 
back  -  Background  Color. pixel ; 
fore  -  Foreground  Color .pixel ; 
font  -  Body  Font  Info; 

max  per  row  -  (  (rlght  edge  -  leftedge) /wva 1 )  -  3; 

1 f ( f%max _per_row  —  0)  j 

no  rows  -  1 /max_per_row  -1; 

|  else  ( 

no  rows  -  1/max  perrow  ; 

)  /*if*/ 


for (1-0; l<-no_rows; !♦*)  ( 

lf(l--no  rows  (4  Umax  per  row  !-  0)  ( 
no  col  -  Umax  per  row; 

)  else  ( 

no  col  -  max  per  row; 

) ✓* 1 f */ 

. . ******* 

Do  a  shell  sort  on  each  row  using  x-coordlnate 


1m  -  1  'max  per  row; 
for (gap-no_col/2?  gap>0;gap/-2) 
for ( }-gap; j<no  col ; j4  4 ) 

for (k-)-gap? (k>-0  44  xleft{k+lm)  > 

xlef t ( k+ lm*gap| ) ;k--gap) ( 
klm  -  k  ♦  lm; 

tempx-xleft (k lm| ; 
tempy-y le  ft | k 1ml ; 

t  emps -aqua  res  f  k  1  m ) ; 

xleft (klm)-xl eft (klm*gap| ; 
yleft[klm)-yleft (klm* gap) ; 
squares (klm) -squares (klm^gap) ; 
xleft [k lmfgapl -tempx; 
yleft ( k lmtgap} -tempy; 
squares [k lm*qapl -temps; 

) /* for  k*/ 


Try  to  eliminate  crossing  lines 

. **»**../ 

/*fp  -  fopen ("debug",  "w")  ;  */ 

1 1  o  h  [  0 1  - 1  ; 
ht  o 1 1 0 1 -0; 
prevx  -  x left ( lm) ; 
no  same  -  0; 

f lag  -  0; 

samexl (no  same]  -  0; 


for (k-0;k<no  col;k++) l 
klm  -  k  +  lm; 

xpolnt  fk  lm  |-CetProfileV|ewX(view,  ProfileView, 
xleft (klm) ) ;  # 

ypolnt (k lm) -Get Prof l leVlewY (view, Prof 1 le  View, 
yleft (klm)); 

If  (klm  <  max_per  row) ( 

yva 1 ( klm)  -  bottom  edge  -  TEXTJKARGIN 
-  Body  Font  In fo->helght ; 
xval(klm)  -  left  edge  4 (wval*wval/3) *klm 

4  wva 1/5  4  wval*  (max  per  row  -  no  col) /2; 

|  else  ( 

yvalfklrol  -  top  edge  ♦  (klm/maxperrow  -1)  * 

(BodyFont Info->helght  4  TEXT  MARGIN); 
xval(klm|  -  left  edge  4  (wva 1 *wva 1 / 3) * 

(klm%max_per_row)  4wval*  (maxper  row  - 
no_col)/2  4  wva 1/5; 

I 

lf(  k ! -0  44  xleft (klm) 'prevx  <  0.001)  ( 
samexr ( nosamej -k lm; 

if  ((1--0  44  xpolnt(klm)<  xva l | klm) ) 1 | 

( 1  > 0  44  xpolnt(klra|>  xval 

1 1 oh ( no  same |  -  1 ; 
htol (no_same|  -  0; 

I 

If  ((1--0  44  xpolnt(klm|>  xva I ( k 1 m ] )  I  1 

(1>0  44  xpolnt ( k Ira) <  xval 

It oh (no  same)  -  0; 
htol (no  same)  -  1; 

I 

flag  -  1; 

1 f (k  —  nocol  —  1)1 
no  same* 4; 

1 

)  else  ( 

If  (flag  --  1)  nosame**; 

flag  -  0; 

ltoh( nosame) -1 ; 

ht  o 1 ( nosame | -0 ; 

samexl ( nosame)  -  klm; 

prevx  -  xl eft ( k lm | ; 

I 

)  /*for  k*/ 

for(n  -  0;  n  <  nosame;  n*4)( 

/ *  f pr lnt  f ( f p, "% 1; xl-%1 ,  xr-% l \n“, n, samexl (n| , samexr In));*/ 
for (gap- (samexr l n ) -samexl (n) 4 1 ) /2;gap>0;gap/-2) 
for ()-gap; )< (samexr (n| -samexl (n)+l);)++) 

f or (k- )-gap; (k>-0  44  ( (yleft (k*samexl (n| * lm)  > 

yleft ( k4 lm*  samexl (n| 4gap) ) 44 ( It oh ( n  |  --1 )  | ) 

(yleft (k4samexl (n)*lm)  <  yleft(k4lm 

4 samexl (n) 4gap) )44(htol(n)--l))  )  ;k--gap) 
klm  -  k4 lm4sam«xl (n ) ; 

/* f print f (fp, "swap  %1  %1;  htol  %l  Itoh  %  1 \n", klm,  gap.  It  oh (n | ,  ht ol [ n ) )  ;  * / 

tempx-xleft (klm) ; 
tempy-yleft (klm); 
temps-squares (klm); 
xleft (klm)-xleft [ klm+gap) ; 
yleft (klm)-yleft ( k lmtgapl  ; 
squares [k lm| -squares ( k lm* gap) ; 
xleft (k lm*gap| -l  empx; 
yleft (k lm*qap) -t empy; 
squares (klm* gap] -t  emps; 
xpo lnt  (klm* gap)  -xpo  1  nt  (klm); 
ypolnt  ( k lm* gap) -ypo! nt  (klm); 
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|  /*for  n‘/ 
/•close (fp) ;*/ 

| /‘for  !•/ 


xpolnt  (klm) -Get Prof  1  leVlewX  (view.  Profile 
xleft ( k lm| ) ; 

ypolnt  (k  Ira )  -Get Prof  1  leViewY  (view.  Profile 
yl  eft  |  k  lm) ) ; 


I 

/ 


for(lndex-0; indexcl; Index**)  ) 

sprint  f (value, "I f", squares (Index) ) ; 


prlntText (vlew,xval [ Index) ,yval ( index) , 
wval, hval, value,  font, fore, 
back,  CLIP  RIGHT) ; 

If (Index  >-  max _per_row  ||  l<-max_per_row)| 

yval( Index)  -  yval( Index)  ♦  Body  Font  In fo->helght; 
) 

XLl ne  (view. self, xval( index) ,yva 1 1 Index) , 

xpo 1 nt ( 1 ndex ) , ypo 1 nt ( 1 ndex ) ,  1 , 1 , 

(back  More)  ,GXxor,  A1 IP  lanes) ; 


XFluah  0  ; 


)/*for  loop*/ 


Create  Racple  Input 

This  routine  reads  a  ploygon  in  the  SIMPL-DIX  linked  list 
and  decomposes  It  Into  an  upper  and  lower  layer.  The  layers 
are  written  to  a  file  which  can  be  used  as  Input  to  RACPLE. 


Revised  10/13/88  EH S 


Bugs:  Sometimes  it  falls  to  print  out  the  second  layer 
I  think  there  is  a  bug  In  one  of  the  possible 
Cases . 


•••*•*•••/ 

Creat  eRacplelnput  (f  1  lecnt ,  Poly,  xlef t ,  zieft) 
1 nt  fl lecnt; 

slmplPolygon  *Poly; 
float  *xleft, *zleft; 

( 

FILE  * fp; 

float  Path  *pathptr; 

f  !  oat  xml  n,  xmax,  zmin,  zmax,  xmaxz,  xmlnz; 
f  1  oat  xma  x zm,  xm  1  n  zm; 
float  xhead,zhead; 
f  loat  x()|  500)  ,  zO f  500)  ; 
f loat  x 1 ( 500  | , zl  (500) ; 
float  x?| 500) , z2 ( 500| ; 

I  nt  1 , 1 , k  ; 

1 nt  head  loc,deqen  edge; 
char  f 1 1 ename ( 80 | ; 

1 nt  count , ( xmln, I xmax, I xmaxm,  Jxml nm; 

1  nt  number  pt  si,  number  pis?; 


/ . 

find  x  and  z  extrema 

*  * . / 


xmln  -  1000.0; 
xmax  -  -10.0; 
xma xz  -  -10.0; 

xmlnz  -  -10.0;  , 

xmaxzm  -  1000. 0; 
xmlnzm  -  1000. 0; 

zmin  -  1000.0; 
zmax  -  -10.0; 
count  -  0; 

pathptr  -  Poly->path; 
xhead  -  pat hptr->polnt .x; 

zhead  -  pat hpt r->polnt . y; 
while  (pathptr  !-  NULL)  ( 
xO(count)  -  pathptr->polnt  .x; 
zO(count)  -  pathptr->polnt . y; 

1 f (pathpt r->point .x  >-  xmax) ( 

If  (count  !-  0  t(  pat hpt r->polnt ,x  —  xmax) | 

1 f (pathpt r->polnt .y  >-  xmaxz) ( 

Ixmax  •  count; 

xmaxz  -  pathptr->polnt  .y; 

I  else  If (pat hpt r->polnt .y  <-  xmaxzm)) 

1 xmaxm  -  count; 

xmaxzm  -  pat hpt r->pol nt .y ? 

I 

)  else  | 

Ixmax  -  count; 

1 xmaxm  -  count; 
xmax  -  pathpt r->polnt .x; 
xmaxz  -  pat hpt r->polnt .y; 
xmaxzm  -  pat hpt r->pol nt .y; 

I 

)/*lf  pathptr  •/ 

1 f (pathpt r->pol nt . x  <-  xmln) ( 

if  (count!-  0  4&  pathptr->polnt .x  —  xmln) ( 

1 f (pat hpt r->poi nt .y  >-  xmlnz) I 
lxmln  -  count; 
xmlnz  -  pathptr->polnt.y; 

)  else  If (pat hpt r->polnt .y  <-  xmlnzm)  ( 
lxminm  -  count; 
xmlnzm  -  pathpt r->pol nt . y; 

I 

)  else  ( 

lxmln  -  count; 

lxminm  -  count; 

xmln  -  pathpt r->polnt .x; 

xmlnz  -  pat hpt r->polnt . y; 

xmlnzm  -  pat hpt r->pol nt .y; 


)/Mf  pathptr*/ 

1 f (pathpt r->point .y  >-  zmax) 

zmax  -  pat hpt r->polnt .y; 
i f (pathpt r->polnt .y  <-  zmin) 

zmin  -  pat. hpt  r->po  1  nt  .y; 
pathptr  *  pat hpt r->next ; 

count ♦  ♦  ; 

)  /‘while*/ 


determine  relative  location  of  head  vertex 
and  cor  ne  r  pol nt  s . 

. . 

1 f ( l xml n-- 1 xm l nm  l xmax-- 1 xmaxm) ( degen  edge  - 

l  else  If  (lxmln  -  -  lxminm)  (degen  edge  «=  1; 

|  else  If  (Ixmax  --  1 xmaxm)  l degen  edge  -  ?; 

)  ♦*  1  se  I  degen  edge  -  <>;  | 


i; 
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If  (ahead  --  xmln  41  (head  —  xmlnz)  ( 
headloc  -  0; 

)  el  se 

If  (ahead  --  xmaa  ft  ahead  —  xaiaxi)  ( 
head_loc  -  2; 

|  else 

If  (ahead  —  ainln  ft  ahead  —  xmlnzm)  I 
headloc  -  i; 

I  else 

lf(xhead  --  xmax  ft  ahead  --  amaaam)  ( 
headloc  -  4; 

|  else 

lf(  lxmln  >  lxmax  ft  laminin  >  lamaam  )( 
head  loc  -  1; 

)  else 

1  f  (  lxmln  <  lamaa  ft  lxmtnm  <  lxmaam  )( 
headloc  -5; 

)  else 

1  f  (  degen  edge  —  0  II  degen  edge  —  1)1 

If (ahead  —  xmax  ft (ahead  <  xmaxa  ft  zhead  >  xmaxzm) )  ( 
headloc  -  3; 

I 

|  else 

1 f (  degen  edge  —  0  I  I  degen  edge  —  21 ( 

lflxhead  --  amln  ft  (zhead  <  xmlnz  ft  zhead  >  xmlnzm) )  ( 

head  loc  -  >; 

I 

I 


switch  (head  loc)  ( 
case  0: 
case  6: 
case  7; 

numberpt  s  1  -  lamax  -  lxmln  il; 
number  pts2  -  lxmlnm  -  lxmaxm  *1; 

1 f (head  loc  —  6|  ( 

number_pts2  ■  count  -  lxmaxm  +1  ; 

I 

for  (k-0;k<number_ptsl;k*‘)  I 

xl  (k)  -  aO |  laml n«k  |  ; 
zl(k|  -  zO(lxmln«k|; 

I 

for (k-0;k<number  pts2;k«*l  I 

1  f  (head  loc— 6)  laml  nm-count  ; 
x2 ( k |  -  a0( lxmlnm-k I ; 
z2 ( k |  -  aO ( lxmlnm-k I ; 

I 

1  f  (head  loc  —  6)  I 

a2(0|-x0|0|; 

z2[0|-z0|0|; 

I 

break ; 
case  2: 
case  3: 
case  4: 

/*  needs  patch  for  case  2‘/ 
number  pt  si  *  lamax  -  lxmln  *1; 
number  pt  s?  *  lxmlnm  -  lxmaxm  *1; 
f or  (k -0;  k<number  pts2;ki'l  I 
x? | k |  -  xO | laml nm-k I ; 

7? \k\  -  z0| Ixml um-k | ; 


for  (k-0;k<number  ptsl;k**)  ( 
xl (k|  -  aoi lxmlntkl ; 
zl | k |  -  zO ( lxmlntkl ; 

I  . 

break ; 
case  1: 

number_pts2  ■  lamtrun  -  lxmaxm  +  1; 
numberptsl  -  lxmax  ♦  count-lamln  ♦  1; 
for  (k-0;k<number_pts2;k+t)  ( 
x2(k)  -  xO (laml nm-k I; 
z2|k|  -  zO[lxmlnm-k | ; 

I 

for (k-0;k<number_ptsl;k+*)  I 
If (k  <  count-lamln) I 

xl ( k ]  -  a0| lxmlntkl ; 
zl ( k ]  -  a0( lxmlntkl ; 

)  else  | 

al|k)  -  x0|k- (count-lxmln) 1 ; 

zl(k|  -  zO (k- (count-lxmln) 1 ; 

I 

I 

break; 
case  5: 

number_ptsl  -  lxmax  -  lxmln  +  1; 

number  pts2  -  lxmlnm  t  (count-lxmaxm)  *  1; 
for  |k-5;k<number  ptsl;k*t 1 ( 
xl | k |  -  xO | lxmln«k| ; 

zl | k |  -  z0| lxmln*k| ; 

1 

for  (k-0;k<number_pts2;kt  t )  ( 

If (k  <-  lxmlnm) | 

x?|k|  -  xO | lxmlnm-k ) ; 
z2 | k |  -  z0( lxmlnm-k | ; 

)  else  ( 

x2(k)  -  xO ( count- (k-lxmlnm) | ; 
z2(k|  -  zO|count- (k-lxmlnm) 1 ; 

1 

I 

break ; 
default : 

number  pt  si  -1 ; 
numberpt  s  2  - 1 ; 
xl (0) -99. 99; 
zl | 0] -99. 99; 
x2 (0| -99. 99; 
z  2 | 0 | -99. 99; 

break ; 

)  /‘switch*/ 

•xleft  -  xmln; 

•zleft  -  xmlnz; 

/•*** . . . 

/*  write  layer  Information  to  fTf  format  file  V 

............ . . . . 

sprint  f  (filename, "RACPLE  Input .11", fl lecnt) ; 
fp  -  f open ( f 1 1 ename, "w" ) ; 

/•••••“write  first  layer  to  f  1  le •»•*•***•***«•* / 

AddPot  ntstol.ayer  (x2,  z?,  inumber  pt  s2)  ; 

fpr  int  f  (fp,  "\n  %10.6f  %10.6f  %10.6f  %10.6f\n", 

xmln, xmax, zmi n- zmax-0 .1,0.11 
fprlnt f (fp, \n  ",2.0>; 
fpr  int  f  (f  p,  "%10.f;(  \n",  1  .0*  number  pts2); 
for  ()-0;  )<.numbi>r  pt  r,7;  )*  '  ) 


I 
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fprintf  (fp,M10.«f  %10.6f\n-,x2(}|, 

z2  (  j )  -zmax  ); 

/*  ••••••write  second  layer  to  file**********4**/ 

AddPolnt  stoLayer  (xl,  zl, (number  _ptsl)  ; 
fprintf)  fp,  M10.6f  \n*,1.0*number_ptsl) ; 
for  (}-0; }<number_pt »1; jss) ( 

f print  f(fp,M10.6f  »10.6f\n*,xl(ll, 

zl ( jl-zmax  ) ; 

I 

fclose (fp)  ; 

I /‘Create  Racple  Input*/ 

AddPolntstoLayer (xx, zz,nn) 
float  *xx,*zz; 
lnt  *nn; 

( 

float  eps; 
float  mlnlength; 
lnt  1, 1,  f  lg.dbgcnt; 
float  xnew( 500| , znew) S00| ; 


This  routine  call*  RACPLE  with  the  appropriate  Input 
file.  Output  la  sent  to  RACPLE_Output . • 

The  Output  file  Is  read  and  the  number  of  lateral  ( 
squares  for  the  first  layer  Is  returned 


revised  10/13/88  EWS 

****** . . 

RunRacple  (fllecnt,output_data) 
float  ‘output  data; 
lnt  fl leent; 

I 

FILE  *fp; 

FILE  *fp2; 
lnt  ldum; 

char  str(17|,str2(7|; 

char  cotrmandl  80 ) ; 
char  Racple[80); 
char  debug) 80] ,  out f 1 le) 801 ; 
char  messagest r lng ( 8 0 1 ; 


mlnlength  -  f abs (xx [ *nn-l I  -  xx[0|)/20.0; 
eps  -  0.0001; 
xnew|0|  -  xx(0|; 
znew(0|  -  zz(0|; 

1  -  i; 

dbgcnt  -  0; 
whlle(dbgcnt<10|| 
dbgcnt  *  *; 

1-1; 

flq  -  0; 

for (1-1; l<*nn; 1  ♦♦) I 

1 f (fabs  (zz( 1 1 -zz| 1-1 1 1 <eps  St  tabs (xx 1 1 1 -xx 1 1 -1 1 1 >mlnlength) ( 

fig  -  i; 

xnew) ) |  —  ( XX ( 1 ) +xx (1-1 )  1  /2.0; 
znew( j|-(zz(l|‘zz(I-l))/2.0; 

1“; 

I 

xnew) j|  -xx ( 1 1 ; 
znew) II -zz 1 1 1 ; 

)**; 

I 

1  f ( f 1 g  —  111 

for (1-1; 1< j; 1 I 

xx | 1 ) -xnew) 1 | ; 

iz( 1 | -znew | l | ; 

/•““patch  to  eliminate  small  sharp  spikes 
from  RACPLE  Input  “““*/ 

1  f (1<  1-11  ( 

1 f  (fabs (xnew) 1* 1 1 -xnew) 1-1 1 )  <  ml nlengt h/ 3 . 5  ts 

( znew ( 1 -1 | <  znew 1 1 | ) St 
(znew (1*1 | <znew I  1  1 1 1  I 
zz(ll  -  ( znew | 1 -1 | * znew ( 1 U | I /2.0; 

I 

I 

1 

*nn  -  1; 

I  else  return(0); 

I 

1 


sprint f (message_strlng,  "Runnl ng  RACPLE  :  %1 * ,  f 1 leent  1 ; 
Prompt (messagest  r 1 ngl ; 


1 f (RACPLE  Path  !-  NULL! ( 

st repy (Racple, RACPLE  Path) ; 

|  else  | 

strepy (Racple, DEF_RACPLE  PATH) ; 

) 

/‘strepy  (Racple,  ‘/users2/edschecli/DEVICESIM/ racple/ racple-)  ;*/ 
sprint f (command, -»s  <  RACPLE  I nput . %l  >  RACPLE  Out  put .% 1  - , 
Racple, f 1 leent ,  f 1 leent )  ; 


system (command) ; 

/•  sprint  f  (debug, "debug.%1",  f 1 leent) ; 
fp2  -  fopen (debug, -w- ) ; 

fprintf  (fp2,  Ms",  command)  ; 
fclose(fp2);  */ 

spr lntf (out f 1 le,  -RACPLE  Output . »1 ", f 1 leent) ; 
fp  -  fopen (out f 1 le,  "r") ; 

f  scanf  (fp.  Ml  »17c  %f  » 7c",  t  ldum,  st  r,  output  data,  st  r7) ; 
/•fscanflfp,  Mf",  out  put  data)  ;  */ 
fclose (fp) ; 

)  /‘RunRacple*/ 


. . * . / 

D1 spl ay  Nonplanar  Et ch ( 1 

. . *““““/ 

I 

dlxVlewport  view; 

Pixel  fore, back; 

Font  Info  ‘font; 

FI  I.F.  *  fp,  •  fopen  ()  ; 

float  xmx, xmn , ymx, ymn, nl ns; 
float  no  pt.  s; 
float  xxl , yyl , xx2,yy?; 
lnt  1,  1.xl,yl,x7,y7; 

vlew  -  DIX  Viewport 1 1  I ; 
bark  -  Background  Color. pixel; 
foie  -  Foreground  Co  I  or . p 1 xe I ; 
font  -  Body  1  ont  I  n  1  o; 


/ 
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fp  -  fopen  (“SAMPLE_netchf27*,  ”r")  ; 
f scant  (fp,  "»f*f%f»f*,  f  xnn,  ixmx,  symn,  tymxl ; 
f scant (fp, Mf, tnlns) ; 
for ( j-0; j<nlns; }♦♦) I 

/*•*** 

read  a  line  segment 

*•***/ 

fscanf  (fp,  *t  f,  tnopts)  ; 
f scant  (fp,  *%Mf*,txxl,  Syyl) ; 
yyl  -  (Prof 1 leBound. top  r  yyl); 

*1  -  Getprofl leVlewX  (view, Prof  1 le  Vlew, xxl ) ; 
yl  -  GetProflleVlewY(vlew,Proflle_Vlew,yyl); 
for  ( 1 -1 ; 1 cnopt s; 1+*) ( 

fscanf  (fp,  Mf%  f“,  txx2,  (yy2)  ; 
yy2  -  (Profile  Bound. top  +  yy2(; 

x2  -  GetProf 1 leVlewX (view,  Prof  lie  View, xx2) ; 
y2  -  Get  Prof  1 leVlewY (view,  Prof  1 le  View, y y 2 ) ; 

/***• 

draw  It  to  screen 

XLlne (vl ew.se  If, xl,yl,x2,y2, 1,1,  (back A fore) ,GXxor,Al 1  Planes) ; 
xriushO  ; 

xl  -  x2; 

yl  -  y2; 

) /*  f or  IV 
) /* for  j V 
f close ( fp) ; 


) 
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/*F77Layer«. c 

*  Routines  used  to  run  non-plsnar 

*  etch  simulation  with  SAMPLE 

* 

‘Edward  H.  Scheckler  Dsc.  1,  1901 

•Copyright  (C(  1988  U.C.  Berkeley  SAMPLE  Group 

•/ 

•include  “SIMPL.h" 

•Include  "SIMPLMacros. h" 

•Include  “Local. h“ 

•Include  <stdlo.h> 

typedef  struct  layers  I 
float  x [ 500|; 
float  z(500|; 
char  layer_type(15| ; 
lnt  number_pts; 

)  LAYERS; 

static  LAYERS  layerarray (30) ; 
static  lnt  numberlayers; 

static  POLYGON? TR  topjpolyl  1st [30|  ; 

/••• . “••••**••***•/ 

Run  nonpl anaretch  0 

. . . . 

I 

FILE  * fp; 

char  coimandllOO) ; 
char  Samplel_7 [ 100| ; 

GLprlnt f ("Initial lrlng  for  nonplanar  etch*); 

1  f  ( -1  I-  Wrlte  SAMPLEnplnrO)  I 
sprint  f (Samplel_7, " /users 2/edscheck/bln/sampl el .  7") ; 

•lfdef  UNIX 

GLprlntf (“running  SAMPLE1.7  .  .  ."); 
sprlntf  (command,  "9s  <  SAMPLE  nonplanar  >  SAMPLE  nOutput ", Samplel  7) ; 
system (comnand) ; 

spr 1 nt f (command, “mv  -f  f77punch7  SAMPLE  net chf 77"); 
system(command) ; 

•endl f 
I 
I 

Run  nonpl ana r  update (flag) 
lnt  flag; 

I 

FILE  • fpt r; 
char  YORN | 1 I ; 
char  conmandl 100) ; 

If  (flag  —  1)1 

Redef  Ine  Proflle  () ; 

I  else  | 

fptr  -  fopen ("tempaave", "r") ; 

Read  Dat  a?  ( fptr)  ; 

(close (fptr); 

I 

spr lnt f (command, " rm  -f  temp  save"); 
syst  em (command) ; 

I 

. . * . . 

Redefine  Profiled 

. . **•**•/ 

1 

r<X)RI)l  NATES  I  p(  1000  I  ,  t>t  (  1  000  )  ; 


lnt  ntp, nbt; 
lnt  1  ,); 

FILE  ‘fp; 

float  xmx, xmn, zmn, zmx, nlns, nntp; 


fp  -  fopen ("SAMPLE_netchf77",  “r"l  ; 
f scant (fp, "9f9f9f9f", i xmx, t xmn, t zmn, t zmn) ; 
f  scant  (fp,  "%f",(nlns); 
for ( j*0; j<nlns; )tt)  I 

fscanf (fp, "9f", tnntp) ; 
ntp  -  l*nntp; 

for (1-0; 1 <ntp; ltt) I 
If  (J— 3)  I 

fscanf (fp, "9f9f",tbt(l| ,x,4bt (1|  .z)  ; 
nbt  -  ntp; 

I  else  ( 

fscanf (fp,"9f%f",4tp|l|.x,ttp|l|.z); 
I 
I 
I 

fclose (fp) ; 

/•  for (1-0; l<layer_array (3 | .numberpts; 1 *<) ( 

bt[l].x  -  layerarray ( 3 | . x ( 1 | ; 
bt[l|.z  -  layerarray ( 3) . z [ 1 | ; 

) 

nbt  -  layer  array(3) .number  pts;*/ 


/*  If  (layerarray  [2|  .layer  _type  —  NULL)*/ 

sprint  f  (layer  array (2 1 . layertype, "PSG") ; 

New  Top  Poly (tp, ntp, bt, nbt, layer  array  I  2) . layer  type) ; 

|/*Redeflne_Proflle*/ 

. . . . *••••••••/ 

GetLayers  () 

This  routine  extracts  the  profile  stored  In  SIMPL-2  and 
Interprets  It  as  a  series  of  layers  which  can  then  be  written 
to  a  SAMPLE  plot  file  In  f77punch7  format. 

It  Is  also  used  to  write  an  Input  file  for  SAMPLE  nonplanar 
etch 


( 


BUGS  -  destroys  SIMPL-2  data  structure 

This  has  been  remedied  by  saving  profile  In  a  data  file 
so  that  It  can  be  reloaded  at  the  end  of  this  routine. 

FILE  *f pt r, • fopen ( | ; 
lnt  1,); 

POI.YGONPTR  top  polygon,  tmpl  pol ygon,  top2  pol ygon; 
POl.YGONPTR  Find  Polygon  ()  .Find  top  polyd; 

VERTEXPTR  top  v,top?  v; 

VERTKXPTR  tmpl  v; 

VERTEXPTR  Get  LeftTopd; 
char  “Write  saveflled; 
cha r  pr  ev I ous ( 1 S | ; 
char  Materia]  name | IS  I; 
chat  *savtt  nattK*; 
cha  r  s  name  I  ft  I)  I  ; 
ctiar  YdRN  |  I  I  f 
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char  layerf Ue_nam«[  14) ; 
float  xleft  (30), xrlght  (30) ,  deft  (30)  ; 

float  topi,  top2  jt;  /*  highest  point  In  layer  */ 
float  bottom  of  poly  0 ; 


save  SIMPL  data  structure  since 
subsequent  code  will  destroy  It. 


fptr  "  fopen("temp_save","ww); 

Write  savefi le2 (fptr) ? 
close (fptr) ; 

The  Layers  are  extracted  by  reading  the  top  contour 
then  deleting  the  polygon  which  make  up  the  top 
contour  which  also  shares  the  same  material  type  as 
most  vertex.  This  continues  until  the  substrate  Is 


the  top 
reached. 


1-0; 

for  (;;) 

I  3-0; 

top2_z-0.0; 

top  i-O.O; 
top  polygon  -  NULL; 
top2  polygon  -  NULL; 
tmpl  polygon  -  NULL; 

tmplv  -  GetLeftTopO ; 
top2_v  -  tmplv; 
t op_v-tmpl_v; 

while  (tmpl _v->xi .x  <  xmax)  ( 

1 f (st  remp (tmpl _v->bMtrl, "AIR") --0)  1 

layer  array |1J ,x[ JJ  -  tmpl_v->xz.x; 
layerarray ( 1 ) . z( j++ 1  -  tmpl_v->xz . z 
strepy (previous, tmpl _v~>aMtrl) ; 

Move (ttmplv, tmpl  v->aMtrl) ; 

I  else  ( 

layerarray(l ) .xlj|  -  tmpl_v->xz.x; 
layer  Array ( 1 ) . z { + J  -  tmpl_v->xz.z 
1  f  (stremp  (tmpl_v->aMt  rl ,  previous)  --0 
st  repy (previous, tmpl_v->bMt  r 
Move (4tmpl_v,  tmpl_v->bMt 1 1 ) ; 

I  else  ( 

strepy (previous, tmpl_v->aMt  r 
Move (4 tmpl _v, tmpl_v->aMtrl) ; 


If (tmpl_v->xz.z  >  top  z)  | 

topz  -  tmpl_v->xz. z; 
topv  -  tmplv; 

top_polygon  -  Find  top  poly (top  v) ; 

I 

tmpl  polygon  -  Flnd_top_j>oly  (tmpl  v) ; 

1  f  (top  polygon  !  -  NULL)) 
lf(  st rcmp(top_polygon->Name, 

tmpl _pol ygon->Name) ! -0) J 
/*  find  top  of  this  polygon  •/ 

If (tmpl  v->xz.x  >  top?  z)  1 

top?  z  -  tmpl  v->xz.x; 
top?  v  -  tmpl  v; 

lop?  polygon  -  Find  top  poly  (top?  v)  ; 

I 


I 


|  /'while4/ 

/*  read  In  last  vertex*/ 


layerarray | I | .x | ) |  -  t.mpl_v->xz  .x; 
layer  array(l| . z l J ♦+ I  -  tmpl_v->xz . z; 


layer  array  ( 1 1 .  number  jrts  *  J  ; 

I  r  Ibot t  omof  poly  (t  oppol yqon)  <  bottomofpoly  (top2polygon) )  | 
top  polygon  -  top2  polygon; 

I 

st repy (layerarrayl l| . layertype, top_polygon->Name) ; 


1 f (st remp (t op  polygon->Name, “SI "I --0)  ( 
xleft [l|-xmln; 
xrlght 1 1 1 -xmax; 
break; 

) 

Flndxmlnandxmax (toppolygon, Ixleft ( 1 | , t xrlght (l|,izleft(l|); 
top_poly_llst ( 1 |  -  top  polygon; 

DeletePolygonltop  polygon) ; 

l**; 

I 

numberlayers  -  1*1; 

/*  fptr  -  fopen (“temp  save",  "r") ; 

Read  Data? (fptr) ; 
f close (fptr) ;  •/ 

I 

........... 

float  bottomofpoly (plyptrl 

POLYGONPTR  plyptr; 

( 

float  bottom; 


VERTEXPTR  tmpv; 


lnt  1  ; 


If  (plyptr  --  NULL)  ( 

ret  urn (-100.0) ; 

I 

bottom  -  10000.0; 

plyptr->Nvertex  -  CountVert 1 ces (plypt r) ; 
tmpv  -  plyptr->HeadVertex; 
for  (1—1;  1<  (plyptr->Nvertex) ;  1  ♦  ♦ )  I 

1 f (tmpv->xz .x  <  bottom)  bottom  -  tmpv->xz.x; 
Move (ttrnpv, plypt r->Name) ; 

I 

ret  urn (bottom) ; 


/'********♦/ 

POLYGONPTR  Find  top  poly (top  v) 
/•**“****/ 

VERTEXPTR  top  v; 

( 

POI.YGONPTR  Find  Polygon  ()  .top  polygon; 
char  Material  name(15|; 


if  (Other (top  v ->aMt r 1 , "AIR") --TRUE) ( 
st  repy  (Mat  er  1  a  1  name, top  V“>aMtrl); 

)  else  If  (Other  (top  V  -  >bMt  r  1  ,  "AIR")  -  -TRUE )  ( 
st.  repy  (Mat  or  I  al  name,  top  v-^bMtrl); 

|  else  1  f  (Ot  ho  r  ( top  v-  >rMt  rl. "AIR*)- -TRUE)  ( 
st  repy  (Mat  or  1  a  1  name, top  v->cMtrl); 

)  e  1  so  { 

SlmKxtl  ("error  finding  top  polygon  !ypo'“); 
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l 

topjpolygon  -  FlndPolygon(top_v,  Mater lalname) ; 
ret  urn (top  polygon! ; 

I 

/•**** . . . . *“*/ 

Write  SAMP LEnplnrO 

. . . . . . . 

< 

lnt  1, j; 

FILE  *fp,  *fp2; 

float  rate,  steps, starttime,etchmodel; 
lnt  lsteps, lstartt lme, letchmodel ; 
char  askstrlng(80|; 

Get  Layers  () ; 

. . . . 

/*  write  SAMPLE  Input  file  for  nonplanar  etch  */ 

/*•**• . . . . 

If  (numberlayers  >  5)  ( 

GLprlntf ("Too  many  layers  for  this  version  of  SAMPLE"); 
return (-1 ) ; 

I 

GLprlntf ("creat lng  nonplanar  Input  file"); 

fp  -  fopen ("SAMPLEnonplanar", *w“) ; 
fprlntf (fp,  “etchnumlay  %d  ;\n",number_layers) ; 

for  ( 1  -number_layers-l ;  1  > - 0 ;  1--I  ( 

1  f  ( 1  —  0  )  ( 

fprlntf (fp,  “etchprofXn") ; 

I  else  | 

fpr lnt f (fp, "nonplanar  Ad\n*,number_layers  -  1  -1  ); 

I 

for  ()-0; )<layer_array( 1 | .  numberpt  s; J**) ( 

1 f (I ( J !-0  is (layer  array |1 I .  x( 11 --layer  ar ray | 1 |.«|)-1| 

ll  layer  array 1 1 1 . 1 |)| --layerarray [1 ) . i ( J-l J ) ) )  ( 
fprlntf (fp, "A10.6f  A10.6f\n", layer  array! 1 1 .x| J| , 

-1 . 0* (layer  array | 1 | . z| j)  -zmax)  ); 

) 

I 

fprlntf  (fp,"  ;  \n"); 

I 

/•  spr lnt f (askst ring, ’Etchmode  1  (see  SAMPLE  user  guide)"); 

Ask  (askstring)  ; 

Answer  Float  Proc (askst  rl ng, let chmode 1 ) ; 
letchmodel  -  1 "etchmodel; V 
letchmodel  -  1; 

f pr 1 nt f ( f p, "elchrates  kd  ".letchmodel); 
lord -number_layers-l; l>-0;  1 — )  ( 

spr  lnt  f  (askst  ring,  "etchrate  for  As,  layer  Ad  (um/sec)", 

layerarray (number  layers-1 -1 ). layer  type, 1 ) ; 

Ask  (ask  st  ring)  ; 

Answer  Float  Proc (askst  rl ng, irate) ; 
fprlntf  (fp,  "  %10.6f  ",  rate) ; 

I 

fprlntf  (fp,"  ;  Vn"|; 

fpr lnt f (fp,  “etchplot  1  0  0;\n“); 

I pr 1 nt f (f p,  "et chwl ndow  A10.€f  ; \n", xmax-xml n) ; 

f pr I nt I (f p, "et chaccur  I  ;\n"); 
spr  lnt  f  (askst  r  lng,  "1  Imest  ep  In  seconds" I ; 

Ask  (ask  st  r  1  nq  I  ; 

Answer  float  Proc(askstrlnq,lstarttime); 
spr  I  nt  Mask  st  r  1  ng,  "number  of  steps"); 

Ask  (ask  st  ring)  ; 

Answer  tlo.it  Proc  (askst  r  I  ng,  1st  eps)  ; 
lsteps  -  I  *  st  eps; 


lstartt lme  -  l'starttlme  ; 

fprlnt f  ( fp,  "etcht lme  Ad  Ad,  Ad; \n", lstartt lme, 
lstartt lme* lsteps,  lsteps) ; 
fprlntf (fp, "etchrun  ;\n"); 
fprlntf (fp, "end  ;\n’l; 
fclose(fp); 

GLprlnt  f ("SAMPLEnonplanar  written") ; 

) 

Find  xml n  and  xmax (polypt r ,  x le  ft ,  xr 1 ght , z left  I 

Routine  to  find  the  leftmost  point  and  right  most 
point  of  a  polygon  and  return  the  x-coordlnate  values 
of  those  vertexes 

struct  Polygon  *poly  ptr; 
float  *xleft, ‘xrlght, *zleft; 

( 

struct  Vertex  ‘tempVertex; 

•xleft  -  10000.0; 

•xrlght  -  0.0; 

•zleft  -  0.0; 
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Catalogue  of  useful  C  functions  in  SIMPL-DIX  and  SIMPL-2  listed  according  to  their  use. 


Catalogue  of  useful  C  functions  in  SIMPL-2 

This  compilation  lists  several  functions  in  SIMPL-2  which  are  useful  for  designing  interfaces  with 
other  programs.  A  brief  summary  of  the  purpose  of  each  function  is  listed  along  with  its  location  in  the 
SIMPL-2  source  code.  The  version  heading  indicates  whether  a  function  is  include  in  the  release  ver¬ 
sion  (R),  or  was  developed  subsequently  by  Alex  Wong  (ASW)  or  by  Edward  Scheckler  (EWS).  For  a 
complete  description  of  the  function  consult  the  comments  in  the  source  listing  for  SIMPL-2.  The  func¬ 
tions  -listed  here  are  grouped  in  the  following  categories: 

Traverse  a  string,  locate  something  in  the  data  base 
'Alter  the  data  base 
Test  for  a  condition 
Utility  functions 

Prompts  and  communication  with  user 
Load/Save 


Function  Name 

Location 

Version 

Traverse,  locate 

GetJLeftTopO 

Returns  a  pointer  to  the  vertex 
at  the  top  left  of  the  profile. 

String.c 

R 

Get_RightTopO 

Returns  a  pointer  to  a  vertex 
at  the  top  right  of  the  profile. 

String.c 

R 

Find_Polygon(V  ertex,Name) 

Returns  a  pointer  to  a  polygon  of  a  given  name 
which  contains  the  specified  vertex. 

String.c 

R 

Move(Vertex,Material) 

Moves  to  the  next  vertex  on  the  polygon  with 
the  given  material  name. (Moves  clockwise). 

String.c 

R 

Move_Backward(V  ertex,Material) 

Moves  counterclockwise  around  the  polygon 
to  next  vertex  with  given  material  name. 

String.c 

R 

Gei_LeftPolyO 

Returns  pointer  to  the  top  left  polygon. 

String2.c 

ASW 

Get_RightPolyO 

Returns  pointer  to  the  top  right  polygon. 

String2.c 

ASW 

Find_Last_PolygonO 

Returns  pointer  too  the  last  polygon  in  linked 
list  data  structure. 

String3.c 

ASW 

Fmd_topjpoly(V  ertex) 

Returns  pointer  to  the  polygon  with  the 
given  vertex  as  long  as  one  of  the  material 
names  stored  in  that  vertex  is  AIR 

F77Layers.c 

EWS 

-  2  - 

Get_LayersO  F77Layers.c  EWS 

Extracts  profile  and  stores  it  in  an  array  as  a  series 
of  layers. 

Find_xmin_and_xmax(Polygon,xl,xr,zl)  F77Layers.c  EWS 

•Finds  x  and  z  coordinate  of  leftmost  vertex 
and  x  coordinate  of  rightmost  vertex 

Alter  Data  Base 

Delete_Polygon(Polygon)  String.c  R 

Removes  polygon  from  data  base. 

Separate_Polygon(Polygon)  String.c  R 

Removes  polygon  pointer  without  changing  the 
linked  vertexes. 

Insert_Polygon(Polygonl  Folygon2)  String.c  R 

Inserts  a  polygon  node  Polygon  1  following 
the  node  Polygon2 

Append_Polygon(Polygon)  String.c  R 

Append  Polygon  to  last  polygon  in  data  structure. 

Polyize.AirO  String3.c  ASW 

Makes  AIR  into  a  polygon 

Rename_Polygon(old_name,name, Vertex)  String3.c  ASW 

Renames  a  polygon 

Merge_PolyO  String3.c  ASW 

Merges  all  polygons  which  touch  and  are  of  same 
material  type. 

Del_Dup_PolygonO  String3.c  ASW 

Deletes  duplicate  polygons. 

One_D_Search_or_Insert(V  ertex.name.posidon)  String.c  R 

Search  for  a  vertex  with  given  x  position 
if  not  found,  insert  a  vertex. 

Two_D_Search_or_Insert (Vertex^  1  ,z  1  ,x2,z2,p,q,name)  String.c  R 

This  function  scans  through  the  edges  between  vertex 
p  and  q  following  ’name’  and  finds  the  crossing  point 
with  the  line  segment  (xl,zl),(x2,z2). 

Insert(Vertexl,Vertex2,Venex3,name)  String.c  R 

Inserts  Vertex  1  between  Vertex2  and  Vertex  3 
following  name. 

Link  (Vertex  1  jiame.V  ertex2)  String.c  R 

link  Vertex  pointer  1  with  ’name’  to  Vertex  pointer  2. 
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Delete_Link(Vertex,name)  String2.c  ASW 

Deletes  a  link  to  a  vertex  of  given  ’name’ 

Delete_LinktoV(Vertexl,Vertex2)  String3.c  ASW 

Deletes  a  link  between  vertex  1  and  2. 

AIRize_String(Vertexl,Vertex2,name)  String.c  R 

assigns  AIR  to  vertices  along  a  suing. 

Write_a_AIR(vertex)  String.c  R 

Fills  vertex  material  names  with  "AIR" 

Delete_Mtrl(name,Vertex)  String.c  R 

Set  pointer  at  a  give  name  in  Vertex  to  NULL 

Add_Mtrl(Vertexl,Vertex2,name)  String2.c  ASW 

Set  pointer  in  Vertex  1  to  point  to  Vertex  2 
with  given  name. 

Substituted  ertex  .name  1  ,name2)  String.c  R 

Change  name  of  pointer  in  Vertex. 

Renamed  ertex  l,Vertex2,name,name)  String.c  R 

Rename  the  string  between  Vertex  1  and  Vertex  2. 

New_Vertex()  String2.c  ASW 

Returns  a  new  vertex  with  null  pointers 

New_PolygonO  String2.c  ASW 

Returns  a  new  polygon  with  null  pointers 

Copy_Vertexd  ertex  1  ,Vertex2)  String2.c  ASW 

Copies  a  vertex. 

Del_BadHeadVertex(array_of_vertices,integer)  String3.c  ASW 

Cleans  up  a  list  of  vertices  with  a  bad  head  vertex. 

Redefine_ProfileO  F77Layers.c  EWS 

Takes  last  contour  in  SAMPLE_netchf77,  and  uses  it 
to  update  the  profile  after  a  nonplanar  etch  simulation. 

Not  yet  fully  implemented. 

Test  for  information 

Count_Vertices  (Polygon)  String.c  R 

Counts  the  number  of  vertices  in  a  polygon 

Count_Polygons(Root_Polygon)  String.c  R 

Counts  the  number  of  polygons  in  the  data  base 

Which_Nodedertexl,Veriex2,name)  String.c  R 

Find  pointer  in  vertex  containing  ’name’ 

Inside_Polygon(flagPolygon,x,z)  String.c  R 
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Determines  if  point  is  inside  given  polygon 


Is_There(Name,  Vertex) 

Checks  if  there  is  ’name’  in  vertex. 

String.c 

R 

Material  (Vertex  .name) 

_ Checks  if  there  is  an  adjacent  material 
which  is  not  AIR  or  BNDR 

String.c 

R 

OtherOtame  1  jiame2) 

Checks  if  material  name  is  other  then 
nil,  BNDR  or  AIR 

String.c 

R 

Verify(Poiygon,name,error) 

Traces  through  the  polygon  to  test  if  the 
polygon  is  well  formed.  Error  string  set 

String.c 

R 

Full(Vertex) 

String.c 

R 

Returns  TRUE  if  vertex  pointers  are  all  assigned 

Thickness  (V  ertex.nam  e) 

Returns  floating  point  thickness  of  a  polygon 

String.c 

R 

Test_Segment(Vertexl,Vertex2Tname) 

String2.c 

ASW 

Test_Connect(V  ertex  1  ,V  ertex2) 

String2.c 

ASW 

Test_Connect_Mtrl(Vertex  1  ,Vertex2jtame) 

String2.c 

ASW 

Two_Same(V  ertex) 

String3.c 

ASW 

Test  if  there  are  two  materials  of  the  same 
name  in  a  vertex. 


Emp  ty  _  Vertex  (V  ertex) 

String3.c 

ASW 

Tests  if  vertex  is  empty 

Test_Cyc!ic(Vertex  1  ,Vertex2,name) 

String3.c 

ASW 

Utility 

Free_Polygon(Polygon) 

String.c 

R 

Free_Vertex(Vertcx) 

String.c 

R 

Free_Mask(Mask) 

String.c 

R 

Free_Block(B  lock) 

String.c 

R 

Prompts 

GLprintf(char_string) 

Graphics_Utils.c 

R 

Prints  a  character  string  to  the  display 

Ask(char_string) 

Graphics_Utils.c 

R 

-  5  - 


Used  to  generate  a  prompt  to  the  user. 


Answer_Float(requesuanswer) 

Get  answer  from  user,  convert  to  floating  point 

Misc.c 

R 

Answer_Up_Proc(request,answer)  Misc.c 

-Get  answer  string  from  user,  convert  to  upper  case 

R 

Answer_Integerrequest,answer) 

Get  answer  from  user,  convert  to  integer 

Misc.c 

R 

Load/Save 

Write_SavefileO 

Saves  data  base  in  a  file 

Write_Data.c 

R 

Read_DataO 

Loads  data  base  from  a  file 

Read_Data.c 

R 

Write_Savefile(fp) 

Save  data  in  file  indicated  by  file  pointer  fp. 

Write_Data.c 

EWS 

Read_Data2(fp) 

Read_Data.c 

EWS 

Load  data  from  file  indeicated  by  file  pointer  fp. 

A  Few  Useful  Routines  and  Global  Variables  in  SIMPL-DIX 

GetYesOrNo(message)  prompt_control.c  R 

Prompts  the  user  to  click  Yes  or  No  with  the  mouse. 

GetProfileViewX(viewport,profile,x)  view_control.c  R 

Used  to  translate  x  coordinate  to  an  integer 
describing  position  in  the  viewport 

GelProfileViewY(viewport,profile,y)  view_control.c  R 

Similar  to  above 

PrintText(variable  list)  graphics_control.c  R 

SIMPL_PolygonRt 

Global  variable  for  start  of  linked  list 
containing  profile  data  in  SIMPL-DIX 


SIMPL-2  source  files  which  were  altered  for  this  project: 
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Deposition. c 

Do_Process.c 

Etching.c 

Read_Data.c 

Write_Data.c 

F77Layers.c  (new  file) 

SIMPL-DIX  source  files  which  were  altered  for  this  project: 

command_control.c 

dix_main.c 

dix_action4.c  (new  file) 


